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Working toward a model independent understanding of cosmic microwave background (CMB) anisotropies 
^^ \ and their significance, we undertake a comprehensive and self-contained study of scalar perturbation theory. 

Initial conditions, evolution, thermal history, matter content, background dynamics, and geometry all play a 
role in determining the anisotropy. By employing analytic techniques to illuminate the numerical results, we 
are able to separate and identify each contribution. We thus bring out the nature of the total Sachs- Wolfe 



;_( ' effect, acoustic oscillations, diffusion damping, Doppler shifts, and reionization, as well as their particular 

-*— ' ' 

C/3 ■ manifestation in a critical, curvature, or cosmological constant dominated universe. By studying the full 

C^ ■ 

angular and spatial content of the resultant anisotropies, we isolate the signature of these effects from the 



dependence on initial conditions. Whereas structure in the Sachs- Wolfe anisotropy depends strongly on 
the underlying power spectra, the acoustic oscillations provide features which are nearly model independent. 
C^ ' This may allow for future determination of the matter content of the universe as well as the adiabatic and/or 

isocurvature nature of the initial fluctuations. 
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It is only with people who know about the useless, 
That there is any point in talking about uses. 
In all the immensity of heaven and earth, 
A man uses no more than is room for his feet. 
Yet if in recognition we were to cut all else away. 
Would it still be useful to man? 



-Chuang-tzu 



I. Introduction 



With the steadily inereasing number of cosmic microwave background (CMB) anisotropy experiments on 
various angular scales (see e.g. [1]), the empirical reconstruction of the process for structure formation in the 
universe will soon enter a new phase. For this task to succeed, the groundwork for understanding anisotropy 
formation must be firmly laid. While numerical studies of specific models abound, this ah initio black box 
approach is not well suited to the reconstruction problem. One must be able to distinguish between the effects 
of initial conditions, evolution, thermal history, matter content, background dynamics, and geometry. With 
the goal of shedding light on the model independent physical mechanisms involved in anisotropy formation, 
we have undertaken a comprehensive and self-contained study of the scalar perturbations which give rise to 
large scale structure in the universe. 

Of the two general classes of scalar perturbations, the isocurvature mode is by far the less well studied. 
A rich structure of anisotropics under the baryon isocurvature scenario is unveiled by generalizing the original 
model proposed by Peebles [2] to arbitrary thermal histories [3,4]. Yet even the familiar adiabatic case holds 
novel features if one steps beyond the standard fio = 1 Harrison-Zel'dovich cold dark matter (CDM) model 
[5]. In this paper, we extend the highly accurate analytic tools developed for this standard CDM model [6] 
to the general case of arbitrary initial conditions, thermal history, and background dynamics. By employing 
these methods to illuminate the numerical results, we examine the physical mechanism behind the evolution 
of isocurvature and adiabatic fluctuations in an JIq = 1 1 open, or cosmological constant dominated universe, 
allowing for possible late or partial reionization. Focusing on the physical interpretation rather than specific 
model dependent results, we explore the possibilities that these as yet undetermined quantities leave open. 

In §11, we discuss the gauge invariant perturbation equations and their general implications. Unlike 
most previous analytic treatments, e.g. [7, 8], we take a multifiuid approach to realistically describe the 
evolution of each component. Superhorizon evolution, analyzed in §111 brings out the differences between 
the isocurvature and adiabatic modes, including the gravitational redshift [9] and curvature effects [10]. 
Further discussion of open universe peculiarities may be found in Appendix A, and commonly used relations 
in Appendix B. As shown in §IV, intermediate scale perturbations in the photon-baryon fluid evolve as an 
oscillator in the potential well created by the total density perturbations. This leads to the characteristic 
oscillatory "Doppler" peak structure in the CMB at recombination for both modes [11]. Photon diffusion 
however erases these acoustic oscillations at small scales [12]. This is especially important in reionized 
scenarios, considered in §V, where last scattering is delayed and the diffusion length grows to be nearly the 
horizon size at last scattering. In this case, degree scale anisotropics can be dominated by the Doppler effect 
from scattering off electrons [13], which at late times are released from Compton drag. 



Putting these results together in §VI, we examine their imphcations for the observable quantities today. 
By analyzing the full matter and temperature transfer functions, we achieve separation of initial and evolu- 
tionary contributions. Robust features in the anisotropy are singled out as potentially useful for extracting 
information about the background cosmology. We conclude in §VII with some general comments on the 
present status of models given their predictions for CMB anisotropics. 

II. The Evolution Equations and Their Interpretation 

In this treatment, we assume scalar fluctuations about a background Friedmann- Robertson- Walker 
metric, 

-\ 2 

ds^ = 



a 



{-drf+-f,jdx'dx^), (1) 



where c — 1, jij is the 3-metric on a space of constant negative curvature, and rj = J{aQ/a)dt is the conformal 
time. Normalized to unity at matter-radiation equality, the scale factor a evolves as a/ a = Ha/ao where the 
overdot denotes a conformal time derivative and 

H^^(^yi±^noH^,-(^yK+\A, (2) 

Va/l + ao \ a / 3 

is the Hubble parameter with Hq = IQQh km s^^ Mpc^^ as its value today. For spatially flat models, 
the curvature parameter K = —Hq{1 — Hq — 51a) goes zero, where the vacuum density is related to the 
cosmological constant by 51a ~ A/3Hq. 

Perturbations around these background quantities may be represented in various ways under gauge 
invariant theory [14,15,16]. Although all are gauge invariant, they reduce to ordinary perturbation quantities 
for different choices of hypersurfacc slicing [17]. This flexibility in the gauge invariant scheme allows us 
to simplify the physical interpretation. For temperature perturbations, we choose shear free Newtonian 
slicing, since on large scales, they are determined by gravitational redshifts from the Newtonian potential. 
Unfortunately, this choice does not clearly bring out the evolution of energy density perturbations, which is 
best studied in the total matter rest frame representation. To avoid confusion, we will only employ photon and 
neutrino temperature perturbations, O = AT-y/Tj and N = AT^/T^, in the Newtonian representation, and 
energy density fluctuations, e.g. A~^ ~ 5p-y/ p~^ and A^ = 6pi,/ pi,, in the total matter rest frame representation. 

A. The Photon and Neutrino Boltzmann Equations 

The full linearized Boltzmann equation for the evolution of the Newtonian photon temperature pertur- 
bations 0(77, X, 7), is given by [8] 

^(e + vi/) = e + vi/ + i^^(e + *) + f /-(© + *) 
= * - $ + f(eo - e + -i^vi + — 7,7jn;^'), 

where v^ is the baryon velocity (c = 1), ^i{— Xi) are the direction cosines of the photon momentum, So is the 
isotropic component of O, and the anisotropic stress perturbation for the photons 11*^ is defined explicitly 
in Appendix A from the quadrupole moment of O. The last term in equation (3) accounts for Compton 
scattering, where f — Xene<JTO./aQ is the differential optical depth, with Xg the ionization fraction, rig the 



electron number density, and ot the Thomson eross section. The gauge invariant metric perturbations are 
^, the Newtonian potential and $, the perturbation to the intrinsic spatial curvature, which are related to 
the density perturbation through a generalized Poisson equation in §IIC. We will commonly refer to both $ 
and \I> as gravitational potentials. 

If the potentials are static and Compton scattering is ineffective, equation (3) implies 0+^ is a conserved 
quantity. This merely represents what we call the ordinary Sachs- Wolfe (SW) effect: a photon experiences 
a fractional redshift of ^P climbing out of a ^l* < potential well. The effective temperature perturbation 
accounting for this shift is therefore O + $. If ^l* changes, the corresponding gravitational redshift of course 
follows suit. Changes in $ also affect the photons through time dilation. Since these effects accumulate 
along the geodesies, we call the combination the integrated Sachs- Wolfe (ISW) effect. The total contribution, 
derived by Sachs and Wolfe [9] , is a combination of SW and ISW effects and completely describes the effect 
of gravitational redshift on the photons. 

In open universes, the 7^ term in equation (3) does not vanish due to the curving of geodesies. Although 
this would seem to complicate matters, its effect on equation (3) is easy to interpret and compute, once we 
decompose the fluctuation into its normal modes. Plane wave perturbations Q = exp(zk • x), appropriate for 
a flat geometry, must be replaced with the eigenfunctions of the Laplacian for an open geometry [18, 19, 20]: 



V'Q = f'Q\, 



-k^Q, 



(4) 



where "|" denotes a eovariant derivative on the 3-space. Since the eigenfunctions are complete for k > \/—K 
one often introduces the auxiliary variable fc^ = fc^ + K. The subtle question of whether 27r/fc or 2Tr/k should 
be considered as the "physical" wavelength of the mode is examined further in Appendix A. 

Since each eigenmode evolves independently in linear theory, it is sufficient to consider temperature 
perturbations to exist in a single fc-mode,* which can be decomposed into angular moments as [10, 21] 



e(r;,x,7)=^e,(r,)G,(x,7). 



(5) 



e=o 



Here the angular functions Gi are defined in Appendix A such that they reduce to Gi = (— i)^exp(jk 
x)P£(k • 7) in the flat space limit, where Pi is an ordinary Legendre polynomial. 

We can now write equation (3) in the standard hierarchy of coupled equations for the ^-modes: 
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(6) 



- iOe, {£ > 2) 



* As usual, the general case can be recovered by summing a power spectrum of these independent fc-modes. 
It should also be noted that all perturbation amplitudes such as O^ have an implicit fc-dependence. However 
when discussing the evolution of a single fc-mode, we drop the index for brevity. After this section, no real 
space perturbation variables are employed. 



where jivl — VbGi = Vb(— fc)^^7i(5|i, and we have made the replacements such as ^(r;,x) = ^(77)Go(x) = 
^(?7)(5(x) here and below. By analogy to equation (6), we can immediately write down the corresponding 
Boltzmann equation for (massless) neutrino temperature perturbations Af(77,x, 7) by making the replace- 
ments Qi -^ N£,T -^ 0, in equation (6). This is sufRcicnt since neutrino decoupling occurs before any scale 
of interest enters the horizon. 

B. CMB Anisotropies 

Although this Newtonian representation of the Boltzmann equation (6) may cause stability problems for 
its numerical solution [22], it serves to bring out the physics of anisotropies quite well. First, scattering tends 
to isotropize the photons in the electron rest frame, leaving anisotropies only in the unscattered fraction: for 
£ > 2, Qi (X exp(— t), whereas 82 oc exp(— 9t/10) due to the angular dependence of Compton scattering. 
Isotropy also requires V^ = 6i = Vt- Even so, the dipole suffers from gravitational infall due to ^, i.e. the 
SW effect. On the other hand, the ISW effect provides a source to the monopole. 

Since the density of free electrons decreases either due to recombination, or if the universe is reionized, 
to the expansion, the CMB eventually ceases to scatter when the optical depth to the present from Compton 
scattering drops to J ° fdi] = 1. Under the standard recombination scenario, this occurs at z^, ~ 1000, 
whereas for reionized models it is delayed until 

n^h^V^^ f 0.05 ^'/' 



if last scattering occurs before curvature or A domination. 

After z*, the photons effectively free stream to form anisotropies. On the last scattering surface, the 
photon distribution may be locally isotropic while still possessing inhomogeneities, i.e. hot and cold spots, 
which will be observed as anisotropies on the sky today. Free streaming transfers fluctuations to high 
niultipolcs, as the £-mode coupling of equation (6) shows. Consequently, in the absence of sources, the 
monopole coUisionlessly damps. For superhorizon scales fcry ^ 1, the photons can only travel a small fraction 
of a wavelength, and thus the fluctuations remain in the monopole. This is reflected in the fe-dependence of 
this ^-modc coupling. If there is subsequent reionization, superhorizon sized fluctuations will consequently 
not damp by isotropization. 

Due to the more rapid deviation of geodesies, a given length scale will correspond to a smaller angle 
in an open universe than a flat one. Thus the only effect of negative spatial curvature in equation (6) is 
to speed the transfer of power to higher multipoles. Its effect is noticeable if the angular scale 9 ^ £^^ is 
less than the ratio of the physical scale to the curvature radius \J~^K jk. One peculiarity arises however. 
Even for the lowest eigenmode, fc = \/ —K or fc — 0, the i?-mode coupling in equation (6) does not vanish. 
Unlike the flat case, this "inflnite wavelength" mode suffers free streaming damping of low order multipoles, 
once the horizon becomes larger than the curvature radius ri\/—K > 1. The physical origin of this effect is 
discussed further in §VIB and Appendix A. 

Finally, let us state some useful relations. As discussed in Appendix A, the total anisotropy is 

2£+l^ V f dk Me ^3,„ ,2 



4.-^^-2^yT^''^ 



^ ^°° '^^ ^'^^ (l + X/fc2)V2fc3|0,|2, 



27r2 ./,^7f k 2£ 



where the ensemble average anisotropy predicted for an experiment with window function Wi is (AT/T)^ = 
J2{'2i+l)WeCi/4:n with Og evaluated at present. Here Mi = {P - K)...{P - Ki^)/{P - KY and reduces 
to unity in the flat space limit. This implies that the contribution to the anisotropy per logarithmic k and 
£ interval is 

We can also sum in I to obtain 

16 + vl/|2_ ^ IGo + *P + £ i^\^A\ (10) 

£=1 

which measures the total power in a single fc-mode. Since fluctuations are merely transferred to high multi- 
poles by free streaming, this quantity is conserved if <& = "f = f = 0, as is evident from equation (3). 

C. Gravitational Potentials 

We of course have to define the gravitational potentials $ and ^ in order to complete the Boltzmann 
equation (6). ft is useful to introduce the following gauge invariant variables: the total density perturbation 
in the matter rest frame pA^ ~ TliiPi^i^ where the sum runs over all the species present and the density 
perturbations are related to the temperature perturbations by 



(11) 



the total matter velocity {p + p)Vt — '^liPi + P-dV-h where p is the pressure; and the anisotropic stress 
pfl — ^jPilli, where contributions come essentially from the radiation quadrupoles 

n^ = ^92, n, = ^N2. (12) 

Employing the Einstein equations, we may now write the generalized Poisson equation as 

AttG / a ' ^ 
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= 49o + 4^^^ 
a k 


A,: 


-4A^o + 4^-; 
a k 



'^=fc^^3^UJ'^"' 



SttG/ a 



'^+*— ^U'""- 



2 



(13) 



As discussed above, scattering suppresses anisotropics such that 92 — 0, and for perturbations larger than 
the horizon scale, N2 <C iVo- Moreover, in the matter dominated regime, the pressure itself is negligible 
p <^ p. The pll anisotropic stress term can thus be ignored as a first approximation, implying $ ~ — vp. 



D. Matter Components 

The baryons evolve under the generahzed baryon continuity and Euler equations 

(14) 



Ab = -k{Vb - V^) + -A^, 



Vb = --T4 + fc* + t{V^ - Vb)/R, 
a 

where R = ipj^/Ap^ is the scale factor normalized to 3/4 ^t photon-baryon equality. Again if a coUisionlcss 
non-rclativistic particle were present, e.g. CDM or compact baryonic objects [23,24], its evolution would be 
obtained by setting f = 0. 

Well inside the horizon, equation (6) and (11) imply that the photons satisfy a separate continuity 
equation A^ = — YafcKy, which reduces the first baryon equation to the familiar form A;, = — fcVJ,. The 
baryon velocity decays due to the expansion and has a source term from infall into gravitational wells. Thus 
the only effect of the decoupled components is through this potential term. 

Early on scattering makes Vb ~ Vy, which shows that the photons and baryons evolve adiabatically 
A^ = YsAf,, regardless of whether the initial conditions are adiabatic or isocurvature. Yet even if the 
universe remains fully ionized to the present, the baryons will eventually decouple from the photons, because 
f/R = ^/sfip-f/pb) goes to zero in the matter dominated limit. Since t cx ilt,, this epoch is independent of 
rib. The Compton drag on an individual baryon does not depend on the total number of baryons. In fact, 
equation (14) and the Poisson (13) equations show that the drag term cx Vb comes to dominate over the 
gravitational infall term cx fc^ at redshifts above z ~ 200(51o^^) ^ ^e • Thus all modes are released from 
Compton drag at the same time, which we take to be 

Zd = 160(r!o/i')^/-'x-2/5, (15) 

defined as the epoch when fluctuations effectively join the growing mode of pressureless linear theory (see 
§VA). 

It is important to realize that the drag and the last scattering redshift are generally not equal. The 
photons decouple from the baryons before the baryons decouple from the photons in the standard recombi- 
nation scenario. Typically the opposite occurs in reionized scenarios resulting in quite different anisotropies 
for the two cases (see §IV and V). 

We now possess all the machinery necessary to describe the evolution of perturbations. Numerical 
solutions, based on Sugiyama & Gouda [25], are presented in the following sections. However, to shed light 
on these solutions, we also apply analytic techniques in the single fluid (§111), photon-baryon fluid (§IV), 
and diffusive limits (§V). 



III. Large Scale Evolution: Single Fluid Approximation 

Since no causal process such as free streaming or diffusion can separate the components, all fluid velocities 
are equal above the horizon. We can thus describe the coupled multi-component system as a single fluid, 
defined by the total matter variables, whose behavior does not even depend on the ionization history. Its 
evolution is determined by combining the equations for the various species, assumed to be either fully 
relativistic or non-relativistic, i.e. equations (6) and (14) with their decoupled variants, 

At - 3w^At = -(l-^){l + w)kVT - 2 fl - ^) ^wW, 



a \ k"^ I \ k^ J a , , 

/ X (16) 

Vt + -Vt = |-^MAt - (1 - 3w)S] +k^-h(l- ^ ) -^n, 
a 3(1 + w)^ 6 \ k^ J \ + w 

where w — p/ p and S = Am — ^/^Ar, with A„i and A^ being the perturbations in the matter and radiation 
energy densities respectively. As we shall see, S can be interpreted as an entropy fluctuation. In the evolution 
equation for Vt, infall due to the potential ^ is countered by the pressure term A^ at small scales. The two 
are balanced at the Jeans scale. In this section, we solve the evolution equations neglecting pressure and 
anisotropic stress as is appropriate for large scales. 

A. Initial Conditions and the General Solution 

The distinction between adiabatic and isocurvature scenarios lies in the entropy term S of equation (16). 
Its evolution is given by the matter continuity equation (14), i.e. S — k{Vr — Kn), where the matter and 
radiation velocities are defined in a manner analogous to Vt [see equation (11)]. Since all components have 
the same velocity, S" is a constant before the mode enters the horizon and, if it is present, must have been 
established at the initial conditions. 

Under this 5* = assumption, we show in Appendix B that the general growing solution of equation (16) 
is 

At = CaUa + CiUi, (17) 

where the C's are fixed by the initial conditions, and we have neglected anisotropic stress. Although 5* is a 
constant above the horizon, we will define C/ = S{0), in anticipation of horizon crossing. As we shall see 
in the next section, the evolutionary factors Ua and Uj take simple asymptotic forms. However to preserve 
generality, we give the complete expressions here: 

1 



U^ 



„. 2 „ 8 „ 16 16 , 

9 9 9 9 



respectively, where keq = {aH)eq/aa = \/2{aafloH^y/^ corresponds to the scale which passes the horizon at 
equality, and we have assumed 11 = 0. The factor D{a) accounts for pressureless growth 

^^5 r^Wa^ (19) 

4 "" 7 {Ha/aor 

where H is obtained from the Hubble parameter by ignoring the radiation 

"'-(^'"'''«-i^'<-h^ ,20) 



We assume curvature and A dynamical contributions arc only important well after equality a ^ I. Curvature 
dominates over matter at a/a^ > r2o/(l — ilg — ^a)i whereas A dominates over matter at a/ao > (flo/^A)^ 
and over curvature at a/aa > [(1 — flo^ ^a)/^a]^^^- Although we will usually only consider A models which 
are flat, these solutions are applicable to the general case. Before either curvature or A domination, D = a; 
afterwards, it goes to a constant. Moreover if a 3> 1, At oc D regardless of scale and initial conditions as 
discussed below. 

Whereas adiabatic scenarios begin with Ci — and finite Ca, isocurvature universes have C/ ^ 0. If 
only baryons, photons and neutrinos are present, as in the case of the baryonic models, S can be broken 
down into 

S = U (a. - ^A.) + (1 - /.) (^A, - ^A,) ^^^^ 

^ (t ^ Jv)o}]^ + Ju'^bvi 

with an obvious generalization to the case of an additional decoupled non-relativistic component. Here 
the neutrino fraction /,y = p,^/{pv + p^-i) is time independent after electron-positron annihilation, implying 
,fu — 0.405 for three massless neutrinos and the standard thermal history. Since Sh-y = A;,— Y4A^ = S(nh/n^), 
and likewise for the neutrinos, 5*6^ and Sbu represent perturbations to the baryon to photon and baryon to 
neutrino number densities respectively. This in fact justifies our use of the term "entropy" fluctuation to 
describe S. Notice that S = then has an obvious interpretation: since the components cannot separate 
above the horizon, the particle number ratios must remain constant. 

Furthermore isocurvature conditions allow no initial curvature perturbations as the name implies. Thus 
the gravitational potentials ^ and $ vanish at the initial epoch, implying Ca — 0. If the perturbations 
were formed in the radiation dominated epoch, they must be placed in the baryons only. In this case, 
Sb-y = Sbu = S initially, which we will hereafter adopt. The generalization to the case where Sb-y ^ Sbv is 
straightforward. Note that any arbitrary mixture of adiabatic and isocurvature modes is also covered by 
equation (17). 

B. The Evolution of the Fluid Components 

Now let us consider the implications and interpretation of the general solution (17). The results for the 
adiabatic mode are extremely simple. When the universe is dominated by radiation (RD), matter (MD), 
curvature (CD) or the cosmological constant (AD), the total density fluctuation takes the form 
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At/Ca = < a MD (22) 

[ D. CD/AD 

Moreover since S ~ 0, the components evolve together Af, = Ac = ^/i^y ~ V4Aiy where Ac is any decoupled 
non-relativistic component {e.g. CDM). The velocity and potential are given by 

-^{keq/k)(l-3K/P)-^a RD 
Vt/Ca = { - ^{keq/k){l - 3i^/fc2)-iai/2 MD 

- (1 - 3K/P)-^D/k CD/AD 



-l{keq/kfil-3K/P)-^ RD 

*/Ca=<( -|(fce,/fc)2(l-3X/fc2)-l MD 

- jikeq/kfil - 3K/k^)-^D/a. CD/AD 



(23) 



Contrast this with the isocurvaturc evolution, 

( ^{k/keg)^{l-3K/P)a^ RD 
At/Ci^I ^ik/kegy{l-3K/k'')a MD (24) 

[ ±{k/keg)Hl - 3K/k^)D. CD/AD 

From the definition of the entropy fluetuation S (see also Appendix B), equation (18) implies that 

n - fa RD 

Ab/C/^i |[a-i + i(fc/fceg)'(l-3i^/fc2)a] MD (25) 

[ I [a-i + i(fc/fceg)'(l - 3K/k^)D] , CD/AD 

and 

(a RD 

A^/Ci^A,/Ci^ll[-l + ^{k/k,gf{l-3K/P)a] MD (26) 

[f[-l + ^(fc/fceg)'(l - 3i^/fc2)i^] , CD/AD 

for the baryon and radiation components. Lastly, the velocity and the potential also have simple asymptotic 
forms, 

- ^{k/keq)a^ RD 

Vt/Ci= ■[ -2^{k/k,g)a^/^ MD 



-^ik/k,g)D/k,g, CD/AD 



-ha RD 



(27) 



*/^^ = ^ - I MD 

iD/a. CD/AD 

Notice that unlike the adiabatic case Vr and ^ have no explicit curvature dependence in this representation. 
Moreover, although these solutions omit radiation pressure and streaming, they are actually valid for the 
matter all the way to the present if horizon crossing occurs after the drag epoch (see Fig. 1). 

Let us try to interpret these results physically. The isocurvaturc condition is satisfied by initially placing 
the fluctuations in the baryons A;, = C'l with A^ ~ 0, so that At = 0. As the universe evolves however, 
the relative significance of the baryon fluctuation Abph/ pT for the total density fluctuation At grows as a. 
To compensate, the photon and neutrino fluctuations grow to be equal and opposite A-y — A,y — —aCj. 
The tight coupling condition A(, = V4A-y implies then that the baryon fluctuation must also decrease 
so that Afc = (1 — 3a/4)C/. The presence of A-^ means that there is a gradient in the photon energy 
density. This gradient gives rise to a dipole Vy as the regions come into causal contact [see equation (6)], 
i.e. Vy ex krjA-^ ex —ka^Cj. The same argument holds for the neutrinos. Constant entropy requires that 
the total fluid move with the photons and neutrinos Vt = V^, and thus infall, produced by the gradient 
in the velocity, yields a total density perturbation At oc —krj(l — 3K/k'^)VT oc A:^(l — 3K/k^)a^Ci [see 
equation (16)]. This is one way of interpreting equation (18) and the fact that the entropy provides a source 
of total density fluctuations in the radiation dominated epoch [26] . 

A similar analysis applies for adiabatic fluctuations, which begin instead with finite potential ^. Infall 
implies Vr oc krj'^, which then yields At oc —kriVr oc — k^ (1 — 3K / k'^)a^'i' , thereby also keeping the potential 
constant. Compared to the adiabatic case, the isocurvaturc scenario predicts total density perturbations 
which are smaller by one factor of a in the radiation dominated epoch as might be expected from cancellation. 

After radiation domination both modes grow in pressureless linear theory At oc D [c.f. equations (22) 
and (24)]. Whereas in the radiation dominated limit, the entropy term S and the gravitational infall 
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term ^ arc comparable in equation (16), the entropy source is thereafter suppressed hy w = p/p, making 
the isocurvature and adiabatic evolutions identical. Since growth is suppressed in open and A dominated 
universes, the potential $ decays which has interesting consequences for anisotropics as we shall now see. 

C. The Total Sachs- Wolfe Effect: Basics 

As noted in §IIB, the SW effect causes the effective perturbation to be + ^ to account for the 
gravitational redshift. Subsequent changes in the potential of course alter the shift, an effect which is 
approximately doubled by the induced time dilation. This is the fSW contribution. To determine the net 
effect however, we must first derive the value of the intrinsic photon fluctuations O. If fc?7 <C 1, the Boltzmann 
equation (6) reduces to the ISW effect 

Go = -$-*. (28) 

Here we have again assumed 11 = 0, which causes a ~ 10% error. For corrections due to 11 see [6]. 

Since the isocurvature initial conditions satisfy ^(0) = = 6o(0), this implies Oo(?y) = '^(v)- The 
effective superhorizon scale temperature perturbation for isocurvature fluctuations is therefore 

Go + * = 2*. (iso) (29) 

Inside potential wells, the ISW effect makes photons underdense so that the gravitational redshift adds to 
the temperature perturbation. This is a direct consequence of the feedback mechanism which generates the 
potentials (see §IIIB). Note however that in a low floh^ model with standard recombination, the potential 
may not reach its full matter dominated value equation (27) by last scattering (see Fig. 2). 

For adiabatic perturbations, the nature of the growing mode Ug [see equation (18)] fixes the initial 
perturbations to be Oo(0) = —^^(0), reflecting the fact that the photons are overdense inside the potential 
well. Although Ug implies the potential is constant in both the matter and radiation dominated epoch, it 
changes to ^(a) = :^\E'(0) through equality. The ISW effect then brings the photon temperature perturbation 
in the matter dominated epoch to 0(?7) = — |'I'(?7) and the effective superhorizon perturbation to (MD) 

[eo + *] = ^*, (adi) (30) 

which is the familiar Sachs- Wolfe result. 

The above results are valid before last scattering at a* . Again since last scattering often occurs before 
full matter domination, one should employ the full form [9o + '^]('r]*) — 0o(O) + 2'i'(ri^,) — ^(0) instead 
of equations (27) and (30). After a*, the photons climb out of the potential wells, leaving the quantity 
[Oq + ^](?7,t) to become anisotropics today. Additional contributions to the anisotropy may arise during 
free streaming again through the ISW effect. Since photon geodesies are radial in the absence of scattering, 
we may use the radial eigenfunctions of the Laplacian to solve for the free streaming behavior. These are 
denoted by Q(x) ~ A'^(x)y/"(0, </>), where — i^^ — k^/K ~ k'^/K + 1, and the radial distance normalized 
to the curvature radius is x = \/^Kri. The radial function X^^x) reduces to je{krj) in the flat limit. For 
superhorizon scales at last scattering, the resultant anisotropics are (see Appendix A), 

1^ = [60 + nri*)XUx -X*) + /V - mtix - X)dv. (31) 
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The right hand side represents the SW and ISW effects respectively. Since the potentials for both the 
adiabatic and isocurvature modes are constant in the matter dominated epoch, the ISW contribution is 
separated into two parts: 

(a) The early ISW effect due to pressure growth suppression after horizon crossing in the radiation dominated 
epoch. 

(b) The late ISW effect due to expansion growth suppression in the A or curvature dominated epoch. 
We shall now consider these effects in more detail. 

D. The Total Sachs-Wolfe Effect: Detailed Structure 

Equation (31) for the total Sachs- Wolfe effect predicts a rich structure of anisotropics for low Vt^ models 
[5]. However, to build intuition for equation (31), let us first consider the famihar adiabatic J^o = 1 model in 
which the ISW term represents only a small correction [6] . A given /c-mode contributes maximum anisotropies 
to the angle that scale subtends on the sky at last scattering. In the k — I plane, the anisotropy will have a 
sharp ridge corresponding to this correlation (see Fig. 3). Here we have plotted (AT/r)^^,, the logarithmic 
contribution to the anisotropy in k and d defined by equation (9).* The pure spherical Besscl functions ii{x) 
show that the series of ridges and filamentary structures are due to the structure of the radial eigenfunctions 
themselves. Notice that the largest fc-modcs project onto the monopole and do not contribute to anisotropies. 

Now let us move onto the more complicated A and open cases. For A models, the ISW term in equation 
(31) yields both early and late type contributions. As we shall see (see §IV), inside the horizon during 
the radiation dominated era, the potential decays due to pressure. This leads to a significant early ISW 
term which is projected onto a somewhat larger angle than the SW effect itself since it originates closer to 
the present. That this is present before A domination is clear from Fig. 4a. Because the early ISW effect 
approaches its maximum of [$ — <I'](?7*) ~ 2^(ry,t) and the adiabatic SW effect is given by approximately 
1/3^(77*), the early ISW ridge is more prominent than the SW ridge in adiabatic models. However for scales 
that enter the horizon during matter domination, the decay in the potential due to radiation is much less 
significant. Thus the height of the second ridge drops significantly at larger scales (see Fig. 4a). 

After A domination a/ao ^ (ilo/^A)^^'^, the potential once again decays. Note however that for typical 
values of Oq > 0.1, this decay begins only comparatively recently leading to late ISW contributions. This 
has three significant consequences. 
(a) The largest fc-modes contribute little to the anisotropy due to the projection effect. Notice that the late 

ISW contribution intersects the 1 = 2 edge of Fig. 4 at a smaller scale than the SW effect. In Fig. 5a, we 

plot the analytic decomposition of contributions to a /c-mode slice corresponding to these large scales. 

The smaller late ISW contribution in fact partially cancels the SW effect. Since the SW contribution 

has not undergone free streaming oscillations at A domination, the two effects contribute coherently and 

cancel due to the decay of the potential. 



* For representation purposes, we chose the initial weights of the /c-modes to correspond to jCyip ex k 
and |C/p oc k^^ for the adiabatic and isocurvature modes respectively. This docs not sacrifice generality 
since one can easily scale the figure to an alternate initial weighting. Note all contour plots of the anisotropy 
represent the numerical results. 
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(b) Since the potential is still decaying at the present, the late ISW effect can boost the low order niultipoles 
for all scales. In Fig. 5b, we plot a smaller mode and show that the late ISW effect is positive definite. 
Recent contributions have not free streamed to the oscillatory regime. The ridge structure of Fig. 4 is 
due to the late ISW effect adding with every other ridge in the SW free streaming oscillation. 

(c) Contributions are spread out over a time comparable to 770 • As we shall see in §VB, this implies cancel- 
lation of the late ISW contribution as the photon travels through many wavelengths of the perturbation 
[27]. Thus late ISW contributions are rapidly damped as the scale decreases leaving only those scales 
that project onto large angles. 

Together these factors imply that if the fc-modes are equally weighted (scale invariant), the result will be a rise 
toward low multipoles from the late ISW contributions [28] . On the other hand if scales that are superhorizon 
sized at late ISW generation are strongly weighted, there is a relative suppression of low multipoles due to 
SW and ISW cancellation. 

Open adiabatic models follow similar physical principles yet still yield significantly different anisotropics. 
Both the late ISW contribution at large scales and the early ISW contribution at small scales contributes 
near the maximum of 2^(77*). On most scales, the combined ISW effect completely dominates over the SW 
contributions (see Fig. 6). However, just as in the A case, the late ISW contributions boost the anisotropy 
in a larger angle than the SW effect for a given fe-mode. Notice where the late ISW ridge intersects i" = 2. 
For the largest mode k = \J —K^ the SW effect consequently can contribute mildly and cancel part of the 
late ISW effect as in the A case. Moreover, this projection effect implies that at these scales, the late ISW 
effect itself is increasingly suppressed with i (see Fig. 5c) . These curvature scale contributions however are 
suppressed by the cut off in the potential near the curvature scale from the Poisson equation (13). This 
has the effect of converting the SW ridge into a peaked structure and curves the contours of the ISW ridge 
away from fc = \/—K. Of course a change in the underlying power spectrum that weights the /c-niodes can 
partially or completely remove this effect. 

At smaller scales, the late ISW effect completely dominates the low order multipoles (see Fig. 5d). 
Finally notice the evolutionary effect of geodesic deviation. Comparing Figs. 4 and 6, we see that the 
fiuctuations are more rapidly carried to high multipoles than in the A case. 

Isocurvature models differ significantly in that the potentials grow until full matter domination. Strong 
early ISW contributions which are qualitatively similar to the SW term will occur directly after recombination 
and continue until full matter domination (see Fig. 2). Thus the projection of scales onto angles will follow 
a continuous sequence which merges the SW and early ISW ridges (see Fig. 7). 

For the A case, the early ISW effect completely dominates that of the late ISW effect. Thus the analytic 
separation shows that the total ISW and SW effects make morphologically similar contributions and the 
boost in low order multipoles is not manifest. Moreover, the two add coherently creating a greater total 
effect unlike the adiabatic case (see Fig. 5a, b). Open isocurvature models behave similarly except that the 
late ISW contributions near its maximum (late ISW ridge) is not negligible. It is thus similar to the adiabatic 
case (c./. Fig. 5d and 8d) except that it does not usually dominate the total anisotropy. 

Notice also that since there is no curvature cutoff in the potentials, the SW and early ISW ridge extends 
all the way to the largest mode k = \J —K. The projection ridge intersects fc = \J —K at i? ~ 10 which is 
the scale the (i7o = 0.1) curvature radius subtends at a distance 770 (see Fig. 8a and Appendix A). This 
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is indicative of the fact that the lowest eigenmodes k -^ \J —K, fc ^ aU contribute to curvature scale 
fluctuations (see §VIB and Appendix A). 

IV. Intermediate Scale Evolution: 
Acoustic Oscillations and Standard Recombination 

As the perturbation enters the horizon, we can no longer view the system as a single fluid. Decoupled 
components such as the neutrinos free stream and change the entropy fluctuation. However, above the 
photon diffusion scale, the photons and baryons are still tightly coupled until last scattering. Since even 
at recombination, the diffusion length is much smaller than the horizon 77*, it is appropriate to combine 
the photon and baryon fluids for study [30,3f ,32]. In this section, we show that photon pressure resists the 
gravitational compression of the photon-baryon fluid, leading to driven acoustic oscillations [6] 

A. The Acoustic or "Doppler" Peaks 

At intermediate scales, neither radiation pressure nor gravity can be ignored. Fortunately, their effects 
can be analytically separated and analyzed. Since photon-baryon tight coupling still holds, it is appropriate 
to expand the Boltzmann equation (6) and the Euler equation (14) for the baryons in the Compton scattering 
time f^^ [33]. To zeroth order, we obtain Qi = Vy — Vb, and 8^ = for £ > 2. There is therefore no Doppler 
shift from a scattering event, since the photons are already isotropic in the baryon rest frame. Again this 
implies that the photon-baryon evolution is adiabatic A;, = V4^7 even when the general evolution is not. 
The first order expansion yields 

eo + T-^eo + k^clQo = F, (32) 



where the photon-baryon sound speed is 
and the forcing function is 



cl = -^—, (33) 

F = -$ ^— $ *. (34) 

1 + R 3 ^ ' 

The techniques established in our recent analysis of adiabatic fluctuations in a flat universe [6] work equally 
well for the general case. Under this formalism, the gravitational driving forces are treated as known external 
potentials in which the photon-baryon fluid oscillates. The right hand side determines the effect of gravity 
through the ISW effect and gravitational infall fc^^P. Notice that the SW effect due to the photon's subsequent 
climb out of the potential well partially counters infall [6] . The ISW term $ also drives the oscillation and is 
important at horizon crossing for modes that cross during radiation domination. We shall see in §IVC that 
ISW contributions after last scattering can also counter or even overwhelm the SW term. 

The left hand side of equation (32) represents an oscillator whose restoring force is due to the photon 
pressure. This homogeneous F ~ equation can be solved by the WKB approximation, 

Oa = {1 + Ry^^^coskr,, 
eb = {l + Ry^/^sinkrs, 
where the sound horizon is 



(35) 



'^^^^,_2 1 / 6 ^^ VT+R+y^R + R:', ^ ^gg^ 

O kf,q y Req 1 + y' Reg 
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with Req = R{rieq). Noticc that if the sound speed is constant, the dispersion relation becomes lu = kcs 
as expected of acoustic oscillations. The solution in the presence of the source F, constructed by Green's 
method, is [6] 

[1 + i?(77)] 1/490(77) = eoiO)coskrsiv) + ^[60(0) + -7?(0)eo(0)]sinfcr,(r7) 

' drj'[l + Rir]')f/'^sm[krsir]) - kr ^i-q' )]F {-q' ) , 





and fc9i — — 3(G)o + ^)- Although the potentials in F can be approximated from their large (§111) and small 
(§IVB) scale solutions, to show the true power of this technique, we instead employ their numerical values 
in Fig. 9. The excellent agreement with the full solution indicates that our technique is limited only by our 
knowledge of the potentials. In almost all models, the potentials can at least be approximated from the 
calculated matter power spectrum at the present and the general principles of their evolution (see e.g. [6]). 

Some common features of these acoustic oscillations, valid for both isocurvature and adiabatic fluctua- 
tions are worthwhile to note. At the start of the oscillation, the amplitude of the monopole increases with R 
{i.e. Q.jjh'^) due to a reduction in the pressure restoring force. Although both the monopole and the dipole 
subsequently decrease, the dipole does so more rapidly due to an additional factor of r^ = Cg oc (1 + i?)^^/^. 
Thus, when last scattering freezes in the adiabatic oscillations, the temperature fluctuations will be dom- 
inated by the monopole. Furthermore, the amplitude of the monopole is itself modulated since inside a 
potential well, the compressional phase of the oscillation is enhanced and the expansion phase suppressed 
if ^ is comparable to Oo (see Fig. 9b). This can also be viewed as a shift in the zero point of the oscilla- 
tions due to gravity. All even peaks for the adiabatic and odd peaks for the isocurvature models suffer this 
suppression. In models were the pressure is relatively low (high flbh^), the expansion phase may be hidden 
entirely in the final anisotropy spectrum [6] . 

Those features which distinguish isocurvature from adiabatic fluctuations are also apparent by inspec- 
tion. For adiabatic initial conditions, the driving potentials are constant until Jeans crossing, at which point 
they decay (see Fig. 9). On the other hand, for the isocurvature scenario, they grow from zero to a maximum 
at Jeans crossing. Thus the forcing function imitates cos krs and sin kr^ in the adiabatic and isocurvature 
scenario respectively and stimulates the corresponding mode of temperature fluctuations. 

B. From the Jeans to Diffusion Scale 

Well below the Jeans scale, the gravitational driving force can be ignored and the photon-baryon fluctu- 
ations behave as simple oscillatory functions, until the breakdown of tight coupling at the photon diffusion 
scale. At this point, photon fluctuations are exponentially damped due to diffusive mixing and rescattering. 
We can account for this by expanding the Boltzmann and Euler equations for the photons and baryons 
respectively to second order in f^^ [34]. This gives the dispersion relation an imaginary part, making the 
general solution 

Go = ^(1 + Ry^^^Viri, k) cos kr, + B{1 + R)-'^/*V{r], k) sin /cr,, (38) 

where the damping factor is 

15(77, fc) = e-('=/'=°)', (39) 
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with the damping scale 

6 7 T {l + Rf 



k-n' = -Jd,- /' ,r:Z^" ■ (40) 



Diffusion thus dissipates these acoustic waves leading to energy input and spectral distortions in the CMB 
[35, 36]. 

The amplitudes of these oscillations, i.e. the constants A and B, are determined by the total effect of 
the gravitational driving force in equation (32). However, a simpler argument suffices for showing its general 
behavior. As shown in §IIIB, isocurvature fluctuations grow like Aj ~ —aC'i until Jeans crossing. Since the 
Jeans crossing time is aj ^ k^q/k, the isocurvature amplitude will be suppressed by k^q/k. On the other 
hand, adiabatic fluctuations which grow as a? will have a (keq/kY suppression factor. This simple argument 
fixes the amplitude up to a factor of order unity. 

We obtain the specific amplitude by solving equation (16) under the constant entropy assumption 5* = 0. 
The latter approximation is not strictly valid since free streaming of the neutrinos will change the entropy 
fluctuation. However, since the amplitude is fixed after Jeans crossing, which is only slightly after horizon 
crossing, it suffices. Under this assumption, the equation can again be solved in the small scale limit. 
Kodama & Sasaki [8] find that for isocurvature fluctuations. 



A = 0, B = -^^Cz, (iso) (41) 



whereas for adiabatic perturbations, 

if fc ^ fceg and krj^ 1. As expected, the isocurvature mode stimulates the siufcrg harmonic, as opposed to 
cos fcrg for the adiabatic mode. 

We can also construct the evolution of density perturbations at small scales. Well inside the horizon, 
A-y ~ 40o by equation (f f). The isocurvature mode solution therefore satisfies (RD/MD) 

A^/C/^-Vefe") (l + i?)-i/4l?(a,fc)sinfcr,. (43) 

The tight coupling limit implies A^ = ^/i^i which requires (RD/MD), 

Afc/C, = f - ^ f ^ ) (1 + Ry'/'^Via, k) sin fcr,. (44) 



!■) {I + R)-^'^V{a,k)s\nkrs. 



This diffusive suppression of the adiabatic component for the baryon fluctuation is known as Silk damp- 
ing [12]. After damping, the baryons are left with the original entropy perturbation C/. Since they are 
surrounded by a homogeneous and isotropic sea of photons, the baryons are unaffected by further photon 
diffusion. From the photon or baryon continuity equations at small scales, we obtain (RD/MD) 



V,/Ci = V^ICi o, ^ fe") (1 + i?)-3/4p(a, k) cos fcr. 



(45) 
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As one would expect, the velocity oscillates 7r/2 out of phase with, and increasingly suppressed compared to, 
the density perturbations. Employing equations (43) and (44), we construct the total density perturbation 
by assuming that free streaming has damped out the neutrino contribution (RD/MD), 



At/C/ = 



3\/6 fee, 



i?"^(l + Rf/'^Via, k) sin kr. 



(46) 



From this equation, we may derive the potential (RD/MD), 



^/Ci 



3 fk 

'4 



3\/6fce, 1 



R-\l + RY''^V{a, k) sinfcr^ 



(47) 



which decays with the expansion since A^ goes to a constant. In Fig. 10, we compare these analytic 
approximations with the numerical results. After damping eliminates the adiabatic oscillations, the evolution 
of perturbations is governed by diffusive processes. 

A similar analysis for adiabatic perturbations shows that diffusion damping completely eliminates small 
scale baryonic fluctuations. Unlike the isocurvature case, unless CDM wells are present to reseed fluctuations 
(see §V), adiabatic models consequently fail to form galaxies. All adiabatic examples shown here, including 
the open ones, are for CDM universes. 

C. Recombination and Free Streaming 

If last scattering occurs before diffusion has damped the acoustic oscillations in a given mode, e.g. in 
the standard recombination models, they will be frozen into the CMB. A generalization of the free streaming 
equation of (31) gives the resulting anisotropics. 



+ f\^-'^)Xi{x-x')dv'- 



(48) 



Here we obtain the diffusion damped fluctuation at last scattering from equation (37) by the replacements 

[6] 

[Oo + mil.) ^ [Oo + *](r;,)I?(77,, fc), 

(49) 
ei(r7,)^ei(77,)I?(77„fc), 



where the damping factor is averaged over the finite duration of last scattering 



(50) 



Since the visibility function f e^'^ goes to a delta function for large r, this definition also coincides with the 
tight coupling limit, equation (39). For analytic approximations of the visibility function see [6,37]. 

As we have shown in [6], equation (48) describes the final anisotropics due to acoustic oscillations to 
high accuracy for any given model. However, for the task of reconstructing the model from observations, 
it is useful to have a simple estimate of equation (48). As we have already seen with the Sachs- Wolfe 
effect, the presence of the radial eigcnfunctions X^ in equation (48) merely represents the projection of 
a spatial scale onto an angular scale on the sky today. The wavenumbers of the peaks in the spectrum 
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will correspond to the modes in which the monopolc reaches an extrema at last scattering: for adiabatic 
fluctuations fc„j = jTiTr/rg (77^ ) , whereas for isocurvature fluctuations k„i = {m — V2)^/^s('7*)i where m is an 
integer > 1. The flrst oscillation is thus at approximately the sound horizon at last scattering rs{ri^). This 
fluctuation is seen as an anisotropy in the multipole im — Kn^ei which corresponds to the angle subtended 
by the scale km at the distance of the last scattering surface. In the small angle approximation 

rg = — sinh[xo - X*] (51) 

which reduces to rg = rjo — rj^ for a flat universe. 

The rapid deviation of geodesies in an open universe causes a given scale to subtend a smaller angle 
on the sky. In older universes (large present horizon 770 from /i or A), the distance to the last scattering 
surface increases, also reducing the angle a given scale subtends. However, this effect tends to cancel with 
the corresponding increase in the sound horizon at last scattering. The imperfect cancellation in the A case 
pushes the peak to somewhat smaller angles. 

In Fig. 11, we plot the angular location (,p ~ T^rg/rs{ri^) from which the series of peaks can be obtained as 
described above. This simple estimate does remarkably well in tracing the higher peaks which are dominated 
by pure acoustic effects. For the first peak, the potentials are still large enough at last scattering to play a 
subsequent role. As mentioned in §IVA, the potential decays after horizon crossing due to radiation growth 
suppression if the universe is not fully matter dominated. For scales that cross the horizon before last 
scattering, this drives the acoustic oscillations. Afterwards it changes the gravitational redshift that the 
photon would otherwise suffer from the SW effect. We call the latter the early ISW effect. Since the sound 
horizon is always smaller than the particle horizon, it contributes to scales between the first acoustic peak 
and scales which cross the horizon during full matter domination. 

For the adiabatic mode, the early ISW effect partially or fully removes the SW redshift, uncovering the 
intrinsic fluctuation 0o(?y*) from the initial conditions and adiabatic growth, as well as a contribution up 
to <&(77*) from the time dilation. This can be seen in the additional contributions to the rms temperature 
fluctuations after the fluctuation crosses the horizon in Fig. 12a. Since the potential has already decayed by 
curvature or A domination, further decay does not significantly affect the radiation. For scales that cross well 
after radiation domination, the late ISW effect is more significant and serves to distinguish between A and 
open models. At intermediate scales, both effects are important. However, since the late ISW contributions 
are cancelled at small scales (see §VB), the total anisotropy is relatively smaller here. 

For isocurvature models, the potential continues to grow outside the horizon in the radiation dominated 
epoch (see Fig. 12b). Therefore, evolution weights the large scales more heavily than the small. In fact, 
the first acoustic oscillation at last scattering is not prominent due both to this enhancement of the large 
scale and the continuing growth of the potential which enhances the second oscillation. This evolutionary 
tilt toward large scale anisotropies can be countered by changing the initial power spectrum. However, the 
corresponding enhancement of small scale matter fluctuations makes reionization likely in this model. In this 
case, anisotropies are destroyed and then regenerated by diffusive effects. It is to this subject, we now turn. 



18 



V. Small Scale Evolution: 
Diffusion Effects and Reionization 

Below the photon diffusion length, even photon-baryon tight coupling breaks down. Since the photons 
diffuse amongst the baryons, the two fluids must be treated separately. Moreover, if the universe is reionized, 
the diffusion length can grow to nearly the horizon size at last scattering. The window for adiabatic oscilla- 
tions closes, and degree scale anisotropics in the CMB will be dominated by diffusive effects. Isocurvature 
baryon models also behave quite differently from adiabatic CDM models with respect to the matter. For the 
decoupled CDM, density perturbations grow regardless of ionization, providing potential wells into which 
the baryons may later fall. Their absence in the baryonic isocurvature case makes the ionization history a 
crucial ingredient for structure formation under this scenario. Consequently although we retain generality 
for CMB anisotropies, we concentrate on the isocurvature model when discussing the effects of reionization 
on the matter. 

A. Matter Evolution: Compton Drag 

1. Partial or Full Ionization 

Due to the lack of Silk damping, baryon isocurvature models typically have high amplitude small scale 
fluctuations which can collapse immediately after standard recombination at z ~ 1000 [2]. It is possible 
that enough energy is released to immediately reionize some fraction Xe of the electrons. This model will 
effectively behave as if recombination did not occur at all. 

Although the tight coupling approximation predicts Vb and A-y go to zero inside the diffusion length, 
its breakdown keeps this from being exactly satisfied (c./. Fig. 10). The single fiuid Jeans argument of §111 
become invalid. Since in the diffusion limit where S ~ A^ 3> A-y, the effect of radiation pressure on Ay 
in equation (16) is exactly canceled by the entropy term. After complete matter domination, the evolution 
equation and its solution therefore becomes identical to the pressureless case, i.e. all modes grow by the 
same factor D{a) [see equation (19)]. 

We can quantify this with the Compton drag argument of §IIC. After Zd = 160(rio^^)~ Xe , the 
baryons are effectively released from photon pressure. Thus, perturbations will grow such that At (a) = 
[D(a)/ D{ad)\Ci for a ^ ad. An excellent empirical approximation (see Fig. 13) to the behavior at interme- 
diate times is given by 

Afe/C/ = t;(a, Od), 



with the interpolation function 



g(Qi,a2) = l+ „/ \ cxp(-a2/ai), (52) 

D{a2) 



where if ai ^ 02, ^(01,02) -^ D{ai)/ D{a2). The velocity Vt is given by the continuity equation (16). 
Notice that this properly accounts for growth in an open and/or A universe. 
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2. Late Ionization 

Now let us consider more complicated thermal histories. Standard recombination may be followed by a 
significant transparent period before reionization at Zi , due to some later round of structure formation. There 
are two effects to consider here: fluctuation behavior in the transparent regime and after reionization. Let us 
begin with the first question. Closely following recombination, the baryons are released from drag essentially 
at rest and thereafter can grow in pressureless linear theory. The joining conditions then imply that 2/5 of 
the perturbation joins the growing mode [34], yielding present fluctuations of ^ ^/5CjD{z ~ Q)/D{z ^ 800) 
where the residual ionization makes the drag epoch z ^ 800 < z*. The evolution is again well described 
by the interpolation function (52) so that Af,(a) — Q{a,at)Ci. By this argument, the effective redshift to 
employ is zt - V58OO ~ 400 - 500. We take here zt ~ 450. 

Now let us consider the effects of reionization at Zi. After Zi, Compton drag again prevents the baryon 
perturbations from growing. Therefore the final perturbations will be Afc(ao) ~ /S.i,{ai)D{aQ)/ D{ad). Joining 
the transparent and ionized solutions, we obtain 

A,M = (^;"'"*) , "<"« (53) 

\Q{ai,at)Q{a,ad), a > Ug 

which is plotted in Fig. 14. Since perturbations do not stop growing immediately after reionization and 
ionization after the drag epoch does not affect the perturbations, we take Qg = min(l.lai, a^). 

B. Photon Evolution: The Doppler and Small Scale Effects 

We now need to examine the evolution of photon temperature perturbations in light of these results 
for the matter. As the diffusion length overtakes the fluctuation, acoustic oscillations in the photons are 
washed away (see Fig. 15). Since the baryon velocity can grow after Zd so that VJ, 3> V^, Doppler shifts 
off moving electrons will regenerate temperature perturbations. Yet since k > kjj, these perturbations will 
be erased as the photons travel across several wavelengths of the perturbation and are rescattered. Unlike 
the acoustic oscillations, photon evolution before last scattering is inconsequential. This is true even in 
late ionization scenarios. The large acoustic fluctuations frozen in at recombination become anisotropics 
as they free stream to the reionization epoch where they are damped as e~'^ by rescattering. Thereafter, 
fluctuations are regenerated by the Doppler effect at last scattering exactly as in the partially ionized case. 
Doppler anisotropics therefore can be completely described by the matter fluctuations at last scattering [38] . 

Moreover, since the optical depth decreases only due to the expansion, last scattering will extend for 
a period of time comparable to 77*. The later last scattering is, the thicker the last scattering surface. 
Cancellation between photons which last scattered off a crest or trough of the matter perturbation will 
severely damp the Doppler effect on scales smaller than the thickness. Together, this implies that the higher 
the ionization fraction Xe, the more severely damped these anisotropics will be (see Fig. 15). One must be 
careful however to avoid overproducing spectral distortions in the CMB due to scattering off hot reionized 
electrons. Under most plausible ionization scenarios, fully ionized open models are ruled out by the low 
Compton-y distortion [4]. 

We can analytically account for these effects by using the weak coupling approximation [38] which treats 
the photons as diffusing across independently evolving baryon perturbations. Moreover, due to the cancel- 
lation of the Doppler effect, ordinarily negligible contributions become significant and must also be included 
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in this forinalism, e.g. the late ISW effect [27] and second order contributions [3,39]. The dominant second 
order correction, caUed the Vishniac effect [40,41] couples Vf, to the spatial dependence of the scattering 
probability, i.e. since n^ — (rie) (1 + Ab), let f -^ t{1 + A^) in equation (3). Note that this effect is not 
present in the numerical calculation. Ignoring curvature and taking the ordinary Fourier transform, we 
obtain the formal solution for the fcth mode of the Boltzmann equation [27],* 

[6 + *] (r,, k, m) - [e + *] (r/d, k, ^)e*fcM(r;.-r,)g-r(^„r,) ^ [Q^^^ ^ Q^^^ ^ Q^](^^ ^^ ^)^ (54) 

with fc/i = k • 7, the optical depth r(77i, 772) = J^^ Tdrj, Qdsw the Doppler and SW contributions, Qjsw the 
late ISW effect and Qy the Vishniac effect. As noted above, scattering rapidly damps out the contributions 
from before the drag epoch as e^"^. Thus the photon temperature perturbation is a function of the matter 
perturbations alone. These source terms are explicitly given by 

»? 



JVd 

Oiswiv, k, m) - /' 2vi/e-(''^'') e''^^("'-'')rfr;', (55) 



Vd 



Oviv, k, m) = r -'J2 -"'^"(^^ ^') Ab(|k - k'l) fe-^^"''") e^'^^^^'-'^^rfr,', 

where recall that the plane wave decomposition is defined such that 7 • Vft(?7,x) — — i/iVf,(?7, /c) exp(«k • x) 
[see equation (5)]. The visibility function f e""^ picks out the epoch of last scattering, and the second order 
nature of the Vishniac term is reflected in the mode coupling sum. 

For scales smaller than the thickness of the last scattering surface, as determined by the width of 
the visibility function re^'^, the Doppler, SW, and Vishniac effects, will be cancelled by oscillation in the 
integrand of equation (55) for all but the perpendicular /i = mode. Analogously ^e^"^ defines a thickness 
of the "gravitational last scattering surface" under which contributions are also cancelled. 

Linear theory flows are irrotational, 7 • v ex /ifc, and gravitational redshifts are absent in the direction 
perpendicular to the oscillation. Both the Doppler effect and the SW effect thus vanish for /i = 0, implying 
severe cancellation. Because cancellation occurs similarly and |^| ^ |Vb| on small scales, the residual Sachs- 
Wolfe effect will be much smaller that the other two effects. By angularly averaging the first of equations 
(54), we obtain the residual effect on the monopole [38], 

60 ^ Ja, ^ ^f, (56) 

which feeds back through equation (54) into the uncanceled /i = mode [41]. This /i = mode is also how 
all effects avoid cancellation. Small scale density perturbations, with oscillations perpendicular to the line 
of sight, can be in bulk motion parallel to the line of sight. The result is the Vishniac effect: a small scale 
temperature anisotropy due to the increase in probability of scattering off an overdense region. The late 
ISW effect is similar. But note that in the A case, the thickness is comparable to 770 implying cancellation 



* Since the Vishniac effect is not linear, we must consider fc-mode coupling. Therefore, in this and the 
following sections where power spectra are employed, we restore the fc-index of the perturbations to avoid 
confusion. 
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occurs up to scales near the present horizon. Thus whereas in the open case one sees a gentle roll off of 
contributions in £, in the A case, anisotropics fall sharply even from the lowest £. 

This method is valid for calculating these secondary fluctuations for either the isocurvature or adiabatic 
scenario, under any ionization history in which last scattering occurs after the end of the drag epoch.* In 
Fig. 10, we show that this approximation (56) compares well with the numerical result roughly between the 
drag epoch and last scattering, as expected. 

Integration of equation (55) determines the present rms temperature perturbations as a function of the 
underlying matter fluctuations P{k) — |At(77o, A:)p. For late last scattering, the integrands in equation 
(55) are wide bell shaped functions. The functions &dsw and Qjsw are therefore approximately Fourier 
transforms whose contribution to the rms can be approximated using Parseval's theorem [39, 41, 38] 

|e + *lLs(%,^)=^7?^ / dx[\GDSw{x)+Gisw{x)\^ + \Gv{x)\^Iv{k)], (57) 
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where x — rj/riQ. The growth is accounted for by the conformal time integrals over 
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where Do = D{irjo) and the time independent mode coupling for the Vishniac effect is [39] 

, ,,, V {kiiof r . f\,, „, (1 - cos2 0)(l-2j/ cos 0)2 P[fc(l + y2_2y cos 0)1/2] P(fcy) 
iv(k) — 7735 3 — / dy I a[cosO)- 



167r2 773 j^ «y_^ V ' (l + y2 _2j/cos0)2 P{k) P{k) ' 

The Vishniac effect peaks strongly to small scales, whereas the first order Doppler and integrated Sachs- Wolfe 
contribution have the same scale dependence reflecting the cancellation process [27]. 

In summary, cancellation occurs because in the diffusive and free streaming limit photons travel through 
many wavelengths of the matter fluctuation source. Cancellation is particularly severe for the Doppler and 
SW effects due to a lack of a perpendicular mode, but is also present for the late ISW and Vishniac effects. 



* For baryonic compact object dominated models, the density of free electrons may be so depleted that 
last scattering occurs before the drag epoch even with maximal ionization. Though the analysis is more 
complicated, it remains true that scattering attempts to equalize Vb and V^. This boosts Vb and suppresses 
Ky at large scales and vice versa at small scales [4]. 
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VI. Matter & Temperature Power Spectra 

The relative amplitudes of the fc-modes which form the power spectrum are often set by an ab initio 
anzatz. Taking a more agnostic approach which allows for future empirical determination of the weights, we 
examine the transfer functions, e.g. P{k) = |Ay(ao,A:)p — |r(fc)C/(fc)p and |T(fc)C^(fc)p for isocurvature 
and adiabatic matter perturbations respectively. For CMB anisotropics, each €-mode evolves differently and 
thus possesses its own transfer function. We present here the full £ — fc space structure of the anisotropy 
transfer function. To illustrate the effects of altering the fc- weighting, we also present a few specific examples 
for the underlying spectrum. Perhaps the simplest possible choice is a random phase pure power law in k 
initially, i.e. \Ca\'^ « fc" and jC/p ex fc" for adiabatic and isocurvature modes respectively. Although this 
may not be realistic near the curvature scale where geometric effects can introduce novel features [42] , these 
toy models do illuminate the general case. 

A. Matter Transfer Function 

1. Adiabatic Models 

For adiabatic models, the matter transfer function is affected by the dynamics and matter content only. 
Since in low f^o models, matter-radiation equality occurs late, the scale at which the transfer function turns 
over due to radiation growth suppression is larger. Furthermore, growth in the matter dominated epoch is 
suppressed due to curvature and/or A. Combining the standard fitting formula for the numerical results [43] 
with our analysis of growth rates, wc may write the total transfer function as 

T(fc) = D{t^^)^^^^^±^^[1 + 3.89g + {lA.lqf + [b.^qf + [Q.nqf]-^'' (59) 

where q = A:/[f7o^^ t;xp(— 2r2f,)] a-nd is valid for f^t <C i^o. Aside from the small Vl[, dependence to account 
for coupling with the photons, q ex k/k^q. For scales that enter before equality, the perturbations grow as a? 
until Jeans crossing at aj ex (keq/k). Thus the transfer function is flat at large scales and goes smoothly to 
A:^^ at small scales. 

The definition of the adiabatic transfer funtion employed here carries information about the growth 
from equality to the present in the form of D{rio). However, in comparing different Qq models, the epoch of 
equality shifts. A more useful choice requires equal potentials $ at the initial epoch (below the curvature 
scale). Specifically, this amounts to considering the quantity T {k) / (Q.qK^)'^ due to the fc^ from the Poisson 
equation (see Fig. 16). 

This has the added benefit that the /c-space SW contribution will also be the same. If the total ISW 
contributions after last scattering are negligible, this normalization of the transfer function is identical to 
a large angle anisotropy normalization for scale invariant spectra. For tilted spectra, one must account 
for changes in the k to i space projection through tiq — ri^. Since n — 1 0,^ ^ 0.9 models approximately 
satisfy these conditions at the largest scales (see §IIID), the relative amplitude of anisotropy normalized 
matter fluctuations on various scales can be read directly from Fig. 16. Matter fluctuations at the 8/i^^Mpc 
{k ~ O.lft-Mpc^^) scale decrease in amplitude for fixed h due to a change in equality rather than A growth 
suppression [44]. 

For open adiabatic models, the situation is more complicated. As we have shown, the late ISW not the 
SW effect dominates the large angle anisotropics. From this we would expect that the anisotropy normalized 
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matter amplitude would decrease relative to Fig. 16. However, this can be countered by the suppression 
of the gravitational potential (SW and ISW) effects from the curvature term in the Poisson equation (13). 
If the underlying spectrum is taken to be |C^P oc k, these effects in fact nearly cancel. However, we have 
reason to believe that the presence of a curvature scale may influence the initial conditions. For example, 
in the specific open inflationary case calculated in [42], the boost from the late fSW effect dominates and 
further suppresses matter fluctuations. 

2. Isocurvature Models 

In contrast to the adiabatic case, the isocurvature matter transfer function exhibits relatively compli- 
cated structure. On scales larger than the Jeans length, the matter gains a fc^ — 3K tail through the feedback 
mechanism (see §HIB) and grow as D(a) after radiation domination. Below this scale, the perturbations 
have damped oscillations around the initial conditions Cj until the end of the drag epoch. Since the Jeans 
scale goes to a constant in the matter dominated epoch, this implies that the transfer function will have a 
significant peak at the maximal Jeans scale. Note that if the universe is not sufficiently matter dominated 
at last scattering, this scale could be less than its absolute maximum. Thus as the ionization fraction de- 
creases, the peak in the transfer function moves to smaller scales (see Fig. 17). Since isocurvature models 
are motivated by the desire to satisfy observational estimates of Hq ~ 0.2, we will concentrate on the effects 
of ionization history rather than matter content. 

We can in fact deduce some of the properties of the oscillatory regime from our simple analysis. In 
§IVB, we have shown that the oscillations decrease as (1 + R)^^^^keq/k until they absolutely disappear for 
scales smaller than the diffusion length k ^ k]j, leaving a constant tail in the transfer function. However, 
after the drag epoch, all scales grow as D{a) so that the flat tail will have an amplitude which is dependent 
on the ionization history. Notice furthermore that the oscillations become less prominent if last scattering 
is delayed, since both the (1 + R)^^^^ and diffusion suppression increases. 

For A models, the change in the growth rate boosts the amplitude of the transfer function. Since neither 
the Jeans scale nor the drag epoch depends on A itself, the shape of the transfer function is the same aside 
from the lack of curvature effects at the largest scales. Analytic fitting formula may be adapted from the 
fully ionized case given by [45] modified to account for the growth rates presented in §VA. 

Since large scale structure measurements indicate that P{k) ex k^^ at intermediate scales 10^^ < k/h < 
lMpc~^, which fall just below the maximal Jeans scale, the isocurvature scenario must have an initial 
spectrum of n ~ —1 at least at these scales [46]. If the initial spectrum is assumed to be a single power law, 
this implies a very steep matter power spectrum at large scales since At oc (fc^ — 3K)Ci. In other words, an 
isocurvature spectrum with index n corresponds approximately to an adiabatic spectrum of n + 4 at large 
scales, e.g. n — —3 yields scale invariance. The steep n = —1 implies large amounts of small scale power 
which allows for the early collapse of structure and early reionization [2] . 
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B. CMB Anisotropies 
1. Large Angles 

As we have shown in §IIID, the total Sachs- Wolfe effect can lead to interesting structure in the anisotropy 
at large scales in an Oo < 1 universe. But how dependent arc the features on the underlying power spectrum? 
In Fig. 18, we show the full contributions to the anisotropy in both I and k as given by equation (9) for 
adiabatic models. Although we have chosen to represent a |Cyip oc k weighting, any initial power spectrum 
can be obtained from scaling by \Ca\'^ /k. The full information of the two dimensional radiation transfer 
function is contained here. Notice that integration in log A; yields the total anisotropy ex [21 + \)Ct and in 
log^ gives the rms temperature fluctuation for a given fc-mode. 

The adiabatic fio = 1 case shown to full scale in the top left panel of Fig. 18, shows the tight k — i 
correlation of the projection from last scattering (see §IIID). The SW effect contributes at large physical 
scales and the acoustic peaks at small physical scales. An expanded view in top right panel shows the break 
between the two effects around logi* = 1.5 and log(fc*Mpc) = —2 for this model. Pivoting the underlying 
power spectrum around this value of k simply emphasizes one effect over the other. 

The situation is more complicated for A and open models. The ISW term contributes to anisotropies 
for intermediate values of k. For low JIq — 0.1 — 0.3 the early ISW effect fills the gap between the SW ridge 
and acoustic peaks of the f^o = 1 model (see §IIID and JVC). The main contribution comes directly after 
horizon crossing for these intermediate fc-values and thus projects onto lower i?-modes than the SW effect. 

At still larger scales, late ISW contributions become important. For A models, they lead to low t 
contributions since most fluctuations have not had time to free stream to high multipoles and those which 
have are cancelled. For intermediate k, the late ISW effect adds in quadrature to the SW effect. Yet for 
the largest fc-modes, the SW effect itself has not free streamed and the late ISW effect will partially cancel 
it. Thus, depending on the /c-weighting of the initial power spectrum, the late ISW term can have different 
effects. In Fig. 19, we plot the anisotropies for single power law weightings. Notice that the boost in the 
low multipoles from A only occurs for intermediate values of the slope. On the other hand, for isocurvature 
models, the A contributions to the late ISW effect are never prominent due to the dominance of the SW and 
early ISW effects. 

For Vto ~ 0.1—0.3 open adiabatic universes, the total ISW effect almost always overwhelms the SW effect. 
There are two exceptions. Below a certain scale, the late ISW effect is thickness cancelled. Moreover these 
scales are often superhorizon sized at radiation domination so that the early ISW effect does not contribute 
either. At the largest scales, the projection carries the late ISW effect onto the unobservable monopole and 
dipole. Thus just as in the A case, the relative weight of SW versus late ISW increases at large scales. Again, 
the SW and late ISW contributions at the largest scales tend to cancel. This is more important for curvature 
as opposed to A late ISW contributions since the horizon size at curvature domination is smaller than that 
at A domination. Indeed for somewhat higher ilg open models (Oq — 0.5 — 0.8) where the SW and late ISW 
contributions are more comparable, cancellation can lead to a suppression of large angle anisotropies [5]. 

On the other hand, for the largest modes the amplitude of late ISW contribution itself decreases with 
(, due to the projection. Yet to have any net effect, the initial power spectrum must rise sharply to large 
scales to counter the k"^ / K Poisson equation suppression. Even the k^^ rise toward large scales in recent 
predictions of an open inflationary model [42], does not overcome this suppression. Thus it is difhcult to 
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obtain a spectrum with falling anisotropics in open universe; in most cases the lowest order multipoles will 
show a rise in the anisotropy d (see Fig. 19). This is often followed by a dip due to the transition between 
late and early ISW domination. 

Open isocurvature models do not suffer Poisson suppression which makes curvature scale peculiarities 
manifest. Anisotropy contributions come from fc's all the way to the curvature scale k — \J ~K or k — 0. 
Notice that this covers an infinite range in log fc, and yet the contributions retain exactly the same £-space 
structure. As we discuss in Appendix A, the radial eigenfunctions Xl{x) have the peculiar property that even 
as the effective wavenumbcr ly — k/\/—K -^ 0, they possess structure on order the curvature scale and are 
exponentially suppressed thereafter. Although the functions are complete, no random phase superposition of 
them will ever produce structure above the curvature scale. As fc ^ 0, all modes contribute at the angle the 
curvature scale subtends when the anisotropy was generated, e.g. at approximately the distance ryo ^ ^ — % 
for the SW and isocurvature early ISW effect. 

For k-scale invariant potential, random phase weighting, the infinite number of decades in log fc as A: ^ 
causes a divergence in the anisotropy, if no cut off is assumed (see Fig. 20a). Moreover, any spectra that 
places even more power on scales k < \J~-K will result in the same final anisotropy. This peculiarity can be 
seen in Fig. 19 for open isocurvature models with n < —3. Note however that "fc-scale invariance" does not 
imply equal power on all physical scales since all low k eigenfunctions have curvature scale power. Physical 
scales above the present horizon do not contribute to anisotropics despite the apparent divergence from low 
k. For the adiabatic case, the suppression of such scales from the Poisson equation prevents this effect from 
becoming manifest for reasonable n. 

This indicates that for open isocurvature scenarios we must alter the power spectrum from fc-scale 
invariance to have enough power at small scales to form galaxies. For spectra that are strongly tilted to 
small scales, anisotropies converge to approximately £(2£ + l)Ci oc P and become independent of n and the 
model. This occurs for n> \ for isocurvature and n > 5 for adiabatic conditions where recall that there is a 
k'^ difference in the correspondence of n to the matter power spectrum. Because fluctuations are dominated 
by the smallest scale fluctuations present, i.e. those at the photon diffusion length fc/j, equation (8) implies 
that Ci is constant in i as required. For an isocurvature scenario with index —1 < ri < 0, which is of interest 
for structure formation, this asymptotic value has not yet been reached and £{21 + 1)C£ ex £ approximately. 
This corresponds to an effective COBE DMR slope of n^s — 2 [5] implying that isocurvature models have 
significantly steeper anisotropy spectrum than the standard CDM model in which UcS — 1, but not as steep 
as one might naively think. In Fig. 20b, we show such an n = — 1 weighting. Notice that bleeding from 
smaller fc-modes than the main k — £ projection ridge is responsible for filling in the low £ anisotropy. 

In summary, we have identified several independent causes of a downturn of anisotropies at low £: 

(a) The Poisson equation curvature cut off. 

(b) SW and late ISW cancellation. 

(c) Eigenfunction curvature cut off. 

The first effect only occurs in open adiabatic models and manifests itself for Q,q < 0.3. The second effect 
is most significant when the SW and late ISW effects are comparable, e.g. open adiabatic models with 
Q,Q ~ 0.5 — 0.8 [5] and comes from scales which are superhorizon sized at the epoch of late ISW generation. 
The last effect applies if the initial spectrum gives significant weight to randomly phased low k contributions 
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and if the contributions are generated early enough to project onto an anisotropy instead of a nionopole 
fluctuation, e.g. open isocurvature models with n < 3. 

Two effects can give an upturn relative to the underlying power spectrum 

(a) Late ISW contributions. 

(b) High fc-mode power bleeding into low £. 

The late ISW effect predicts a rise toward low £ because of crest-trough cancellation at small scales. In a 
A universe, this cutoff scale is on order the present horizon so that contributions arc already falling sharply 
with £ at low £. For open universes, the late ISW effect contributes earlier and has a smaller scale cutoff. 
Thus the signature of the anisotropy transfer function is a rise to a plateau at low £. However, since the 
/c-modes which contribute to this effect are the intermediate ones, this effect is only manifest if the initial 
spectrum gives them weight. For pure power laws, this requires a roughly scale invariant potential: fc'^$^ = 
constant. In an open universe, the Poisson cutoff can change the plateau to a dip in the anisotropy at low 
£. For the opposite case of small scale weighted power spectra, n > —O.I for isocurvature and n > 3 for 
adiabatic, the higher fe-ridges in the projection effect contribute strongly to low i multipoles. This implies 
that there is a maximum slope with which low order multipoles can rise, Ci ~ constant. 

If any such features arc detected in the observed spectrum and are statistically significant considering 
cosmic variance, some variation of the standard CDM picture will be necessary. However, even though 
a simple tilt (single power law) in the power spectrum cannot mimic such features, it is clear that more 
complicated initial spectra can. This degeneracy between the initial conditions and the evolutionary effects 
is lifted by assuming an ab initio model. In this case, large scale anisotropics are a simple yet powerful probe 
of the underlying cosmology as is well known. Alternatively, once the fundamental cosmological parameters, 
e.g. r^o, h, A, are known, they will tell us what the initial conditions for structure formation are. 

2. Intermediate to Small Scale Anisotropies 

In standard recombination scenarios, acoustic oscillations determine the structure of anisotropies for 
both adiabatic and isocurvature modes. Since these oscillations contain a great deal of structure, it is obvious 
that more cosmological and model information can be extracted here than at larger angles. Moreover, once 
coverage of the sky at these angles becomes more complete, these measurements will be more immune to 
uncertainties from cosmic variance. 

The angular scale of the peaks is determined by the projection of the sound horizon at last scattering 
onto the sky today and is independent of the underlying power spectrum. Three cosmological quantities, 
with their corresponding dependence on fundamental parameters, enter into its construction: 

(a) rs(?7*) ~ f {rioh^ , ilbh^) , the sound horizon at last scattering, 

(b) rjo — rj^, ^ f{iloh^, h, HaH^), the distance to the last scattering surface, 

(c) K — f{Qoh^, h, V,Ah^), the curvature. 

The first task is to distinguish between adiabatic and isocurvature scenarios. For adiabatic models, the peak 
£ values follow the series (1:2:3: 4...), whereas for isocurvature models (1:3:5: 7...). Since the first 
peak is contaminated by the early ISW effect (see Fig. 21), the higher peaks are the most reliable measure 
of this effect. On the other hand, this rise to the first peak can also be used to separate isocurvature from 
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adiabatic models. Wc have noted in §IVC that the first isocurvature oscillation is low in amplitude. Only in 
adiabatic models does the first oscillation truly stand out as a peak. 

Once adiabatic and isocurvature models are distinguished, the location of the peaks is uniquely predicted 
by the cosmological parameters. However, the degeneracy in the dependence on floh^, ^a^^, h, and fife/i^ 
does not allow inversion of the relation [47]. For example, an fib = flo ~ I, h ~ 0.5 adiabatic model predicts 
£ ~ 400 and h — l.O, £ ^ 500 which can mimic projection effects from curvature and A. Of course, if one 
is willing to restrict flbh^ to lie within the nucleosynthesis bounds, its effect on rs{ri^) is negligible. On the 
other hand, the geodesic deviation due to K with low ilg and A == is a severe and easily tested effect. If the 
angular location of the peaks turn out to be multiples of a high £ > 400 — 500, then curvature must almost 
certainly be present in the model since no reasonable change in r,, (ry* ) or r^o ^ V* can account for it [48] . 

Since isocurvature acoustic oscillations are likely to be erased by reionization, let us concentrate on 
lifting the degeneracy for the more plausible adiabatic case. We can use the deviation of the first peak 
from the acoustic series predicted above for this purpose. The early ISW effect pushes the peak to larger 
scales for low ftoh"^ universes. Moreover, the amplitudes of the peaks contain a large amount of cosmological 
information as well. Even though this depends on the underlying power spectrum, a minimal assumption, 
such as a pure power law only over the range of the peaks, would be sufficient to allow interesting constraints 
on cosmological parameters. As we have seen, 

(a) Lowering HqH^ boosts the first peak relative to the higher peaks due to early ISW contributions 

(b) Raising flbh^ boosts the odd numbered peaks over the even due to reduction in the pressure relative to 
the gravitational force. 

In fact, these opposing h dependences nearly cancel for the first peak if ilg = 1 and flbh^ is given by big bang 
nucleosynthesis. This is not true for the higher peaks [6,32]. Thus the relative amplitudes of the series of 
peaks contain crucial cosmological information. These important tests will depend on having experimental 
information for anisotropics £ > 200. 

Considering the present experimental focus on ^ < 200 anisotropics, it would be useful to extract 
information from the ratio of large to intermediate angle anisotropics. For instance, in the n = 1 model of 
Fig. 21a, the rise to the first peak is more dramatic in low flgh^ universes. Unfortunately, this of course 
depends on the specific model in question. However in general, lowering floh^ increases the intermediate 
anisotropics through the early ISW effect whereas increasing flbh^ does the same through the acoustic 
oscillations. For large scales, the late ISW effect can boost anisotropics a comparable amount in the open 
but not the A case. However one must recall that in the open case there is also Poisson suppression of the 
power spectrum and other curvature effects. 

Of course, allowing the thermal history to deviate from the standard recombination scenario introduces 
another degree of freedom which complicates the extraction of cosmological information. If reionization is 
low, the acoustic peaks which are damped as e""^ below the horizon at last scattering, may still be observable. 
However if the ionization is high, the detailed information in the acoustic oscillations is lost to us. This is 
likely to be the case for isocurvature models. If the initial spectrum is chosen to be consistent with large 
scale structure n ~ — 1, the large fluctuations at small scales could result in reionization. Normalized to 
large scale anisotropics, standard recombination models also produce excessively large intermediate scale 
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adiabatic oscillations in the standard recombination scenario (see Fig. 21b). Reionization is therefore also 
necessary. 

In this case the sole feature is the damping scale which measures the photon diffusion length at last 
scattering. In Fig. 21b, we show the effects of altering the ionization history of open and A isocurvature 
models. Assuming a cosmological model, the damping scale fixes the ionization history. On the other hand, 
assuming an ionization history {e.g. fully ionized), it essentially probes the horizon size at last scattering 
as projected via geodesic deviation. Although A models are older and yield a larger distance to the last 
scattering surface, the geodesic deviation effect pushes the damping scale of open models to even smaller 
angles. Notice that this also makes the open universe large angle anisotropics nearly independent of ionization 
history since these angles correspond to supcrhorizon scales at last scattering. 

As for the amplitude of the regenerated fluctuations, we may employ the analysis of §VB, to gain insight 
into the numerical results. In Fig. 22, we show a comparison of isocurvature temperature power spectra from 
the numerical and analytical calculations. The numerical calculations are purely first order and do not include 
the Vishniac contribution. The Doppler and SW (DSW) fluctuations are increasingly suppressed by thickness 
cancellation as last scattering is delayed, as reflected in the time integrals of equation (58) . The late ISW 
effect is of course independent of ionization but increases as ilo decreases. For the fully ionized, low ilg = 0.1 
universe shown here, the late ISW contribution thus more than doubles the temperature fluctuations at 
intermediate scales. On the other hand, the Vishniac effect depends quadratically on the amplitude of the 
matter fluctuations and thus is larger for later last scattering. In our detailed numerical study [4], we show 
how these various effects can be combined to yield the minimal anisotropics for the isocurvature model. 

Reionized adiabatic models look similar to isocurvature models in that the sole feature is at the diffusion 
scale at last scattering. If no underlying power spectrum is assumed, it may be difficult to distinguish between 
the two. However, as large scale structure measurements reach to larger scales and CMB experiments to 
smaller scales, it will be possible to entirely remove the ambiguity of the initial power spectrum (see e.g. 
[49] ) . Consistency between the matter and radiation power spectrum is indeed the ultimate test of any model 
for structure formation. As we have seen, the difference in the matter and temperature transfer functions on 
the same scale can remove all doubt on the question of adiabatic vs. isocurvature initial conditions and/or 
standard recombination vs. reionized thermal histories. 

VI. Discussion 

We have comprehensively studied the evolution of density and temperature perturbations with an ar- 
bitrary spectrum of adiabatic and isocurvature perturbations in a critical f^o — Ij open, and A dominated 
expanding universe. By employing an analytic treatment, we provide model independent insight into the 
formation of anisotropics that is confirmed by its agreement with the full numerical calculation. It thus 
becomes possible to separate and interpret each physical process that generates these perturbations. 

Our treatment identifies numerous sources of anisotropies. Curvature effects due to geodesic deviation 
and on the fluctuations themselves give rise to peculiarities in the anisotropy spectrum which may soon be 
constrained by the observations. Moreover gravitational redshift effects due to the photon's climb out of the 
potential well (SW effect) as well as decay or growth in the potential due to radiation (early ISW effect) and 
the decay due to the rapid expansion in an open or A dominated universe (late ISW effect) carry specific 
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signatures that may be identifiable in the large angle anisotropies. However, the manifestation of these 
effects in a particular model will depend on the initial power spectrum. In examining the dependence on 
initial conditions, we also present a particularly simple derivation of the , (adiabatic) and 2 (isocurvature) 
coefficients multiplying the gravitational potential in the SW effect. 

Smaller angle anisotropies carry information which is less dependent on the power spectrum. We have 
investigated the nature of acoustic oscillations which give rise to peaks in the anisotropy as well as diffusion 
damping which is responsible for its small angle cutoff. Moreover, we have provided a very simple formula 
which predicts the angular location of the peaks as a function of the matter content and geometry of the 
universe. The physical origin of their relative heights is also clarified. In reionized models however, acoustic 
oscillations are damped and give way to last scattering effects due to baryons in infall. At intermediate 
scales, this leads to the Doppler effect whereas at small scales significant second order Vishniac contributions 
must be considered. 

Although the principles outlined here are valid for any model, they can also be used to evaluate currently 
popular models for structure formation. At the present however, it is not even clear which model, if any, is 
consistent with the large scale structure data alone, much less the detailed features in the CMB anisotropies. 
Despite the success of the elegantly simple standard CDM model for structure formation, it is becoming 
increasingly clear that some modification cither in the model or our understanding of its implications is 
necessary (e.g. see [50] for a review). Normalized to large scale anisotropies, standard CDM predicts matter 
fluctuations which imply a moderately anti-biased picture of galaxy formation [51] and more small scale 
power than is observed for peculiar velocities. It is also difficult to understand the dynamical measurements 
of a low ilo at small scales in this picture [52]. The obvious solutions within the context of CDM are to 
either change the initial power spectrum from Harrison-Zel'dovich n = 1, or lower Qq to move the equality 
cut off to larger scales. Indeed the shape of the matter power spectrum alone seems to indicate il.oh ~ 0.25 
[43], and determinations of a high Hubble constant h ~ 0.7 — 0.8, if confirmed, also support low Oq models 
due to the age problem [53]. 

We have fully examined the consequences for anisotropies of these standard solutions. The signature 
of low ilo niodels at large scales depends on the underlying power spectrum. Particularly in the case of 
open models, where we expect deviations from a single power law spectrum at the curvature scale, this 
ambiguity can change the relative amplitudes of anisotropies to matter fluctuations as well as the shape of 
the large scale anisotropies themselves. For A models, this is perhaps less of a concern. The boost in low 
order multipoles from the late ISW effect can be used to constrain n — 1 models [54]. The acoustic peaks 
provide a better handle on the underlying cosmology from both their angular location and relative heights. 
Even with complications such as gravitational wave contributions, which can boost the large scale anisotropy 
relative to the matter [55], the information contained in the acoustic peaks is not lost. 

Another possible alternative is to abandon adiabatic fluctuations in favor of isocurvature ones. This 
model also changes the relative amplitude of matter versus temperature perturbations. However given the 
likelihood of rcionization, the thermal history of baryonic isocurvature models can be adjusted to match 
the observations. The fundamental probe here is the slope of the matter and temperature power spectra. 
Present indications are that n ~ — 1 (ncfr = 2) from large scale structure measurements. The implied 
discrepancy with flat large scale anisotropies with ricff ~ 1 [56] is beginning to indicate that no single power 

30 



law model is adequate [26]. While this is not necessarily surprising for the open version, it would require a 
dramatic break in the power spectrum to counter the heavily small scale weighted power required by large 
scale structure. Perhaps more damaging to this model is the growing body of intermediate scale ^ ~ 50 ~ 200 
anisotropy measurements. If a steep rise toward £ ~ 200 is also confirmed [57], there will also have to be an 
additional break below the curvature scale. Furthermore, there are indications that even large scale structure 
measurements themselves do not fit with single initial power law isocurvature models due to features in the 
matter transfer function [43]. 

Finally a change in the matter content, e.g. adding massive neutrinos [58] or topological defects [59], 
is another possibility. Although we do not explicitly consider such exotic models, the principles outlined 
here remain valid. Sachs- Wolfe contributions and acoustic oscillations are determined from the gravitational 
potential in the same way in these models. Thus once the evolution of the matter is understood, the 
implications for anisotropics is apparent. 

Given that none of these alternatives provide a compelling ab initio model for structure formation, it 
is perhaps best to keep an open mind to all of these possibilities. As the large scale structure and CMB 
anisotropy data continue to accumulate, the general principles formulated here will aid in the empirical 
reconstruction of a consistent model for structure formation. 

Many are those under heaven who attend to their theories and techniques, 
and they all believe that nothing can be added to the ones they possess. 
Where is the true way of old to be found? 

In a haze! Where am I going? 
In a daze! Where shall I arrive? 
With the myriad things before me, 
None will serve as final destination. 

-Chuang-tzu 
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Appendix A: Open Universe Normal Modes 

1. The Radial Representation 

Fluctuations in an open universe must be decomposed in the eigenfunctions of the Laplacian 7*-'Quj = 
—k^Q. To gain intuition about these functions, let us examine an explicit representation. In radial coordi- 
nates the 3-nietric becomes 



7ydxMxJ' = -K-^[dx^ + sinh^ x{dO^ + sin^ Od(P^)], 



(A-1) 



where recall x = V^^^V- Curvature makes the surface area of a shell at distance 77 increase a.s —K ^e^^ 
rather than if' for super-curvature distances x ^ 1. The Laplacian can now be written as 



l'''Q\ij = -i^sinh X 



Tx l^^"^ ^ % 



" ' smt^-— - 



902 



(A-2) 



Since the angular part is independent of curvature, we may separate variables such that Q = Xl{x)Y™{0, 4>) 
where jy^ ~ k"^ /{—K) = —{k'^/K + 1). From equation (A-2), it is obvious that the spherically symmetric 
£ = function is 



x'Ax) 



sin(j^x) 



^-K-^ 



sin(fcA?7) 



(A-3) 



i^sinhx ' ksmh{Ar]^/-K) 

As expected, the change in the area element from a flat to curved geometry causes ^/—Krj -^ sinh x in the 
denominator. The higher modes are explicitly given by [18, 19] 



and becomes jei^kAr)) in the flat space limit. Here 

Mi[k) = - 



(A-4) 



(A-5) 



(fc2 - kY 

which reduces to unity as X ^ 0. It represents our convention for the normalization of the open universe 
functions. 

It is often more convenient to generate these functions from their recursion relations [20]. One such 
recursion relation is 



dr] " 2^+1 



kX, 



e + 1 



i-i{e + 2)^ 



kX, 



(A-6) 



2£+l 

which is of the same form as the Boltzmann equation (6) for {Q + '^)/{2l+l) in the free streaming limit. This 
is quite natural since free streaming photons arrive at the observer on radial geodesies as an examination of 
equation (A~l) shows. Thus the solution of the free streaming Boltzmann equation in the absence of the 
ISW term is obvious: 

9,(7?, fc) 



[eo + *](r/*,fc)x^(x-x* 



(A-7) 



2£-hl 

where we have assumed that the boundary condition at last scattering is given by the monopole fluctuation 
as is appropriate to the SW effect. The ISW effect acts like an impulse (^ — ^)5rj at some intermediate time 
rj which then free streams to the present. The full solution therefore is 



^^ ^ [Go + mv*.~k)K{x - X*) + /'[* - *](r;', ~k)Xi{x - xVv' ■ 



(A- 
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Let us now examine the peculiar nature of the eigenfunctions. Since they are complete for k > \/—K, 
i.e. k > 0, should 27r/fc or 27r/fc be considered the effective wavelength? In Fig. 23, we plot the spherically 
symmetric £ — mode given by equation (A-3). The argument in favor of k is that its first zero is at 
Aj] — Tr/k. This is related to the completeness property: the zero crossing property shows that as A: — > 
we can obtain arbitrarily large structures. However even in this limit, the amplitude of the structure above 
the curvature scale is suppressed as e~^. The effective scale of the prominent structure thus goes to the 
curvature scale favoring k^^ = l/\/—K as the effective wavelength. In fact, the e^-^ behavior is independent 
of the wavcnumbcr and £, if x S> 1. 

This peculiarity in the eigenmodcs has significant consequences. Any random phase superposition of the 
eigenmodes X^ will have exponentially suppressed structure larger than the curvature radius. Even though 
completeness tells us that arbitrarily large structure can be built out of the X^ functions, it cannot be done 
without correlating the modes. This is even if the structure has support only to a finite radius which is 
above the curvature scale. 

Is the random phase hypothesis and the lack of structure above the curvature scale reasonable? The 
fundamental difference between open and flat universes is that the volume increases exponentially with the 
radial coordinate above the curvature scale V{xc) ^ [sinh(2xc) — 2xc] as the line element of equation (A-1) 
shows. Structure above the curvature scale implies correlations over vast volumes [29]. It is in fact difficult 
to conceive of a model where correlations do not die exponentially above the curvature radius. The random 
phase hypothesis has been proven to be valid for adiabatic inflationary perturbations [42]. However, a 
definitive answer to this question for isocurvature models awaits the invention of a mechanism for generating 
such perturbations in a consistent model for structure formation. 

2. General Angular Functions 

Although the radial representation suffices for many purposes, often one needs the full machinery of the 
general normal mode decomposition. Formally, the angular and spatial fiuctuations of the full radiation field 
is decomposed into [10] 

oo 

e(r7,x,7)-^e,(r/,fc)G,(x,7), (A-9) 

i=o 

where 

G,(x,7) = (-fc)-^Q|,^..,,(x)P;-*^(x,7), (A-10) 

and 

^0 = f , Pl^ l\ 



1 

-(07T-7-J, (A_^^) 



n' = 7;{ilW-f'). 






with parentheses denoting symmctrization about the indices. For flat space, this becomes Gg = {—iY cxp(zk- 
x)P^(k • 7), where Pg is an ordinary Legendre polynomial. Notice that along a path defined by fixed 7, the 
flat Gi becomes ji{kri) after averaging over fc-directions. Travelling on a fixed direction away from a point 
is the same as following a radial path outwards. Thus fiuctuations along this path can be decomposed in 
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the radial eigcnfunction. We shall see that this argument can be generalized to the open universe case and 
allows one to interpret equation (A-10) more easily. 

We can also use the properties of Gi to simplify the Boltzmann equation (3). The anisotropic stress 
perturbation of the photons, defined as 



n!/ ^ 4| ^ (^f 7^ - ^7^^) e(ry,x,7), 



is therefore related to the quadrupole moment 

1 



The recursion relation 



,,1^1.^1 = Y^e.G.. 



d 



(A-12) 



(A-13) 



d 



YG.u = -G[x(r7),7(r7)] -X* — G,+f — G, 



drj 
k 



dx^ 



e 



21+1 



\-{e- 1) 



^7* 



(A-14) 



which follows from equation (A-10) and (A-11) [21], completes the simplification of equation (3) to (6). 
Here we take x(?7) to be the integral path along 7. By comparing equations (A-6) and (A-14), the open 
universe generalization of the relation between Gi and the radial eigenfunction is now apparent: 

Gt\An)nij])\=MiXi{ji). 

The only conceptual difference is that for the radial path that we decompose fluctuations on, 7 is not 
constant. This also clarifies the interpretation of the recursion relation for Gg [equation (A-14)]. Finally by 
employing these definitions, we may write the temperature correlation function as [10] 



(e*(77o,x, 7)6(770, X, 7')) = 






'-P\ee{vo,~k)\'Peh-i'), 



(A-15) 



where Pi is a Legendre polynomial. This implies the definition of Cg in equation (8). 

Appendix B. Single Fluid and Other Useful Relations 

Above the horizon the entropy perturbation S is constant, and all perturbation quantities can be 
obtained from the solution for the total density perturbation At. Combining the total continuity and 
Euler equations in (16) yields the second order evolution equation 

2 



d^ 



LA. 

a da 



1 



1 






A7 



k 



1 



■iK 



jS, 



(B-1) 



where 



/ = 



3a 



5 a 



5 = 2 



4 + 3a 2 1 + a 

9a a 6 + 7a 



4 + 3a 2(l + a)2^ 
^2 



(B-2) 



h^- 



3(4 + 3a)(l + a)' 
a 



.? = o 



3(4 + 3a)(l + a)2' 
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where recall that a is normalized to unity at matter-radiation equality. Here we have taken the anisotropic 
stress n = and assumed that the universe is in the matter or radiation dominated epoch. The solutions to 
the homogeneous equation with 5* = are given by 



Ua^ 



,2,8 16 16 , 

9 9 99 



1 



Vr, ' 



"'" + '> (B-3) 



a^/a + 1 



and represent the growing and decaying mode of adiabatic perturbations respectively. Using Green's method, 
the particular solution in the presence of a constant entropy fluctuation S becomes A^ = CaUa + CdUd + 
SUi, where Ui is given by 

4 f^Vf, 3X\ 3a2 + 22a + 24 + 4(4 + 3a)(l + a)i/2 ^ 



^' 15 UeJ V fcV (l + a)(3a + 4)[l + (l + a)i/2]4 "■ (^ ^^ 

After radiation becomes negligible, the both isocurvature and adiabatic modes evolve in the same manner 

/ \ 2 



At + -At = ^nGp — At- (B-5) 

a \aoJ 

For pressureless perturbations, each mass shell evolves as a separate homogeneous universe. Since a density 
perturbation can be viewed as merely a different choice of the initial time surface, the evolution of the 
fractional shift in the scale factor, i.e. the Hubble parameter H, must coincide with At- It is simple to check 
that the Friedman equations do indeed imply 



a ,v . „ / a 



2 



H+-H = inGp { — ] H, (B-6) 

a \aoJ 

so that one solution, the decaying mode, of equation (B-5) is At oc H [34]. The growing mode At oc D can 
easily be determined by writing its form as Z? ex HG yielding 

G+(^+2f)G = (B-7) 



which can be immediately solved as [34] 



^(«) °^ ^ / (^- (^- 



Note that we ignore pressure contributions in H [c.f. equation (20)]. If the cosmological constant A = 0, 
this integral can be performed analytically 



3 ^ 3(l + x)V2 ^^^^^ ^ ^^^^^ ^ ^^^^^ 



D(a) oc 1 + - + "^^^^;/ ln[(l + x)!/^ - x^l^\ (B-9) 



where x ~ (f^Q — l)(a/ao). In the more general case, a numerical solution to this integral must be employed. 
Since before curvature or A domination D oz a, the full solution for At, where the universe is allowed to 
pass through radiation, matter and curvature or A domination, can be simply obtained from equation (B-3) 
and (B-4), by replacing a with D normalized so that D ^ a early on. 
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With the solution for Ay and the definition of S [equation (21)], all component perturbations can be 
written in terms of Ay. For example, in the baryonic isocurvaturc scenario. 



A. = 



1 



4 + 3a 



[45' + 3(l + a)AT] 



and 



A^ = -{Ab-Sbu), 



(B-10) 



(B~ll) 



A^=.-(A,,-^,,^). 



The fact that in this model the curvature perturbation vanishes initially when the universe is radiation 
dominated allows us to set S*;,,/ = S^^. The velocity and potentials can be written as 



3a( 3Ky'l + a_ 
ka\ k^ J 4 + 3a 



dAi 



1 



da 



1 + a 



3 /fceg 

AX k 



1- 



3K 



1 + a 



(B~12) 



A. 



where note that constant entropy assumption requires that all the velocities Vi ^ Vt- The relation for the 
velocity may be simplified by noting that 



rjia)^^[VT+7i-l] 



1 



'-K 



: cosh 



1 



RD/MD 
2(1- Jlo) a 



ilo o-o 



(B-13) 



MD/CD 



where CD denotes curvature domination with A = 0. For A 7^ 0, it must be evaluated by numerical 
integration. Before curvature or A domination 



d (1 + a) 



1/2 



V2a 



(B-14) 



which can be used to explicitly evaluate (B-12). Finally, in Tab. 1 we list some commonly used symbols in 
the paper and the equation in which they first appeared. 



Table 1. Commonly used symbols. Time variables a, z, rj, and x ^re often used inter- 
changably with special epochs listed here under scale factor a entries. Component density 
Ai and velocity Vi are defined in §IIC and D, with i as b for baryons, 7 for photons, i> for 
neutrinos, and c for coUisionless cold dark matter. Note that V^ = 0i (see following page). 
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Symbol 


Definition 


Equation 


At 


Total density fluctuation 


(11) 


e 


CMB temperature fluctuation 


(3) 


Oo 


CMB monopole fluctuation 


(6) 


e, 


CMB i^th multipole fluctuation 


(6) 


n 


Anisotropic stress perturbation 


(12) 


* 


Gravitational (Newtonian) potential 


(13) 


$ 


Gravitational (curvature) potential 


(13) 


* 


Gravitational (curvature) potential 


(13) 


V 


Conformal time 


(1) 


V 


Curvature normalized wavcnumber 


(31) 


(JT 


Thomson cross section 


(3) 


T 


Thomson optical depth 


(3) 


X 


Curvature normalized distance 


(31) 


V 


Diffusion damping factor 


(39) 


G 


Drag growth factor 


(52) 


Ca 


Initial adiabatic spectrum 


(17) 


Ci 


Initial isocurvature spectrum 


(17) 


Ci 


Anisotropy power spectrum 


(8) 


D 


Pressureless growth factor 


(19) 


F 


Gravitational driving force 


(34) 


H 


Hubble parameter 


(2) 


Ne 


Neutrino Cth multipole 


(6) 


K 


Curvature 


(1) 


R 


Normalized scale factor ipb/Ap-^ 


(14) 


S 


Entropy fluctuation 


(16) 


T 


Matter transfer function 


(59) 


Vt 


Total velocity amplitude 


(11) 


xi 


Radial cigenfunction 


(A--3) 


a 


Scale factor 


(1) 


flo 


Present scale factor 


(1) 


ad 


Compton drag epoch 


(15) 


aeq 


Equality scale factor 


(1) 


ai 


Ionization epoch 


(53) 


a* 


Recombination conformal time 


(7) 


Cs 


Photon-baryon sound speed 


(33) 


k 


Laplacian wavenumber 


(4) 


k 


Renormalized wavenumber 


(4) 


ko 


Diffusion damping wavenumber 


(40) 


keq 


Equality horizon wavenumber 


(18) 


i 


Multipole number 


(5) 


re 


Projection factor 


(51) 


Ts 


Sound horizon 


(36) 


Xe 


Electron ionization fraction 


(3) 
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Figure Captions: 

Figure 1. Large scale open isocurvature evolution {Qq — 0.2, h ~ 0.5, no recombina- 
tion). Perturbations, which originate in the baryons, are transferred to the radiation as 
the universe becomes more matter dominated to avoid a significant curvature perturbation. 
Nonetheless, radiation fluctuations create total density fluctuations from feedback. These 
adiabatic fluctuations in A^ dominate over the original entropy perturbation near horizon 
crossing an in the matter dominated epoch. The single fluid approximation cannot extend 
after last scattering for the photons a.^, since free streaming will damp A-y away. After 
curvature domination the total density is prevented from growing and thus leads to decay 
in the gravitational potential ^. 

Figure 2. The total Sachs- Wolfe effect (i7o — 0.1, /i — 0.5, standard recombination). In 
the adiabatic case, temperature fluctuations are enhanced in gravitational wells such that 
and ^ cancel, yielding Gq -I- ^ = Vs^ i^ the matter dominated epoch. For the isocurvature 
case, the ISW effect creates a net total of Gq + ^ = 2^ reflecting the anticorrclatcd nature 
of radiation and total density fluctuations. After last scattering at a^, this SW contribution 
(analytic only) coUisionlessly damps from the monopole. The rms temperature fluctuations 
(numerical only) acquires contributions after a^ from the ISW effect due to the radiation 
(early) and curvature or A (late) contributions. These contributions are relatively more 
important for adiabatic models due to the partial cancellation of Go and ^P at last scattering. 
Since A domination can only have occurred comparatively recently, the late ISW effect is 
also less important in a A compared to an open universe. 

Figure 3. VLq — \ adiabatic full photon spectrum (standard recombination). Shown here 
and in Figs. 4,6,7 is the contribution to the anisotropy per logarithmic k and i interval 
{/yr/T)1^. [equation (9)] with equally spaced contours up to a cut off set to best display the 
features in question. The strong correlation between £ and k merely reflects the projection of 
a scale on the last scattering surface to an angle on the sky. At log^ > 2, SW contributions 
fall off and are replaced by the acoustic peaks (saturated here) . The detailed structure can 
be traced to the radial eigenfunction A^(x) — ji(x) which governs the projection and free 
streaming oscillations. 

Figure 4. A adiabatic photon spectrum (JIq = 0.1, /i = 0.5, standard recombination). 
Unlike the CIq = 1 case, this scenario has significant contributions from after last scattering 
through the early and late ISW effect, (a) The early ISW effect arises if horizon crossing 
is near radiation domination, and projects onto a second ridge which is more prominent 
than the SW ridge at intermediate but not large angles, (b) After A domination, the late 
ISW contributions come free streaming in from the monopole yielding a boost in the low 
order multipolcs for a small range in fc, due to cancellation with SW contributions at the 
largest scales and crest-trough cancellation at smaller scales. Scales depicted in Fig. 5a, b 
arc marked here by dashed lines. 

Figure 5. Analytic separation of adiabatic large angle anisotropics (JIq — 0-1 ^^ — 0-5, 
standard recombination, arbitrary normalization). Scales are chosen to match the features 
in Fig. 4 and 6. A models: (a) At the largest scales, e.g. here k = 10~^Mpc~"'^, the SW 
effect dominates over the late ISW effect due to projection. However since the potential 
decays, the late ISW effect partially cancels the SW effect if the mode is superhorizon sized 
at A domination, (b) Intermediate scale peaks in Fig. 4 are due to the late ISW boost of 
the higher SW free streaming ridges. Open models: (c) The maximum scale corresponds to 
the curvature radius k = \/—K. For this scale, the SW effect projects broadly in i peaking 
near € ~ 10. For the late ISW effect, this scale projects onto the monopole and dipole near 
curvature domination thus leaving the ISW contributions to decrease smoothly with i. (d) 
At smaller scales, corresponding to the large ridge in Fig. 6, the late ISW effect projects 
onto £ ~ 2 — 10 and completely dominates leading to a rising spectrum of anisotropics. 
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Figure 6. Open adiabatic photon spectrum {Hq = 0.1, h = 0.5, standard recombination). 
(a) Like the A case, the radiation ISW effect contributes significantly to intermediate angle 
anisotropics, (b) However, as already noted in Fig. 2, the late ISW effect appearing at the 
left is much more significant than the corresponding A effect. Thus on all angular scales, the 
total ISW contribution dominates the SW effect. Contours curve away from the curvature 
scale log(A:*Mpc) — —3.8 due to suppression of the potentials from the Poisson equation. 
Scales depicted in Fig. 5c, d are marked here with dashed lines. 

Figure 7. Open and A isocurvature photon spectrum {Qq — 0.l,h = 0.5, standard recombi- 
nation). Unlike their adiabatic counterparts, the potential grows in the radiation domination 
era only to turn over and decay in the curvature and A dominated era. The ISW contribu- 
tion will thus smoothly match onto the SW contribution. This has the effect of merging the 
SW and ISW ridges to make a wide feature that contributes broadly in £. For A models, the 
radiation ISW effect completely dominates over the A ISW effect. Scales depicted in Fig. 8 
are marked here in dahsed lines. 

Figure 8. Analytic separation of isocurvature large angle anisotropics {fto = O.I, h = 0.5, 
standard recombination, arbitrary normalization). Scales are chosen to match the features 
in Fig. 7. In general, isocurvature models have strong early ISW contributions which mimic 
and coherently boost the SW effect. A models: (a) Notice that the shape of the SW 
and ISW effects are identical at large scales, (b) Even at the late ISW peak, the early 
ISW contributions are so strong that the late contributions are never apparent unlike the 
adiabatic model. Open models: (c) As with A models, radiation epoch contributions are 
significant making the SW and ISW contributions similar for large scales, (d) Near the 
peak of the curvature ISW contribution however, the relative contributions are similar to 
the adiabatic case. 

Figure 9. The acoustic oscillations (fto — 0.2, /i = 0.5, no recombination). The photon- 
baryon fluid acts like an oscillator in a potential well. The dipole, i.e. the photon velocity 
Vj, is increasingly suppressed with respect to the monopole as (I -t- R)^^'^, where the vS 
accounts for the three degrees of freedom in the dipole. Scales which reach an extrema 
in the monopole at last scattering will correspond to the so-called "Dopplcr peaks" in the 
anisotropy spectrum. Also displayed here is the semianalytic approximation described in 
the text, which is essentially exact. The small difference in the numerical amplitudes of $ 
and ^ is due to the anisotropic stress of the neutrinos. Whereas the isocurvature case has 
Oq = r^b, the adiabatic model has fif, = 0.06 and a consequently smaller R. 

Figure 10. Small scale isocurvature evolution and photon diffusion (JIq = 0.2, h — 0.5, 
no recombination). At small scales gravity may be ignored, yielding pure adiabatic oscil- 
lations. Perturbations in the photons damp once the diffusion length grows larger than 
the wavelength ko < k. Likewise the adiabatic component of the baryon fluctuations also 
damps leaving them with the original entropy perturbation. After diffusion, the photons 
and baryons behave as separate fluids, allowing the baryons to grow once Compton drag 
becomes negligible o > a^. Photon fluctuations are then regenerated by the Dopplcr effect 
as they diffuse across infalling baryons. The analytic approach for the photons in this limit 
apply between the drag epoch and last scattering a^ < a < a^. 

Figure 11. Angular scale of the "Dopplcr peaks" (standard recombination). The physical 
scale of the peaks is simply related to the sound horizon at last scattering. Peaks in the 
anisotropy today will correspond to multiples of the angle that this scale subtends on the 
sky £p = TTrg/rs{ri^), as discussed in the text. Varying A or ft, increases both the sound 
horizon at 77* and the present horizon 770 leaving little effect. For open models, a given scale 
will correspond to a smaller angle by geodesic deviation. This simple analytic estimate for 
the peak location is valid for pure acoustic contributions and underestimates the scale of 
the first peak in low iloh^ models due to neglect of the early ISW effect. 

Figure 12. The total ISW effect {Hq = 0.1, h = 0.5, standard recombination, k = k^ x 
I0~'^Mpc~^). (a) Adiabatic models. The decay of the potential as the scale enters the 
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horizon due to pressure growth suppression causes the early ISW effect which boosts scales 
approaching the first acoustic oscillation. The largest scales which enter after radiation 
domination are boosted by the late ISW effects due to the rapid expansion in open and 
A models, leaving a deficit at intermediate scales, (b) Isocurvature models. Scales which 
enter early during radiation domination do not grow as much due to the suppression in the 
potential. This enhances the large scale with respect to the small. Notice that the second 
adiabatic oscillation (fca = 20) can be comparable to the first since the turnover in $ occurs 
later. Only at the largest scales is the distinction between open and A models manifest in 
the rms temperature fluctuations. 

Figure 13. Open isocurvature baryon evolution under partial ionization (ilo = 0.2, h ~ 
0.5). The baryons are released to grow in pressureless linear theory after Compton drag 
becomes negligible. Since this epoch becomes earlier as the ionization fraction is decreased, 
present day fluctuations will be larger for low Xe , if normalized to the ionization independent 
fluctuations at large scales. Unlike the CDM case, baryons have no potential wells into which 
they might fall after the drag epoch and the transfer function is extremely sensitive to the 
ionization history. 

Figure 14. Open isocurvature baryon evolution in late reionization scenario (fio ~ 0.2, h ~ 
0.5). Here the universe is suddenly reionized to Xg = 1 at redshift Zi after a transparent 
period 1000 > z > Zi. Perturbations are released from drag following recombination only to 
suffer its effects once again between the ionization and drag epochs Zi > z > zj^. Thus the 
final fluctuations will be larger for later reionization. 

Figure 15. Small scale isocurvature temperature evolution under partial reionization (ilo = 
0.2, h ~ 0.5, numerical). If the universe stays transparent after standard recombination at 
z* ~ 1000, the acoustic oscillations in the photon fluid will be frozen. However these large 
fluctuations are suppressed by diffusion damping in partially reionized models. Although the 
Doppler and other diffusive effects regenerate fluctuations at small scales, these effects are 
also suppressed under the diffusion length (i.e. the thickness of the last scattering surface). 

Figure 16. The adiabatic matter transfer function {h — 0.5, fib = 0.01). The transfer 
function has been scaled by {il.oh'^)^'^ to compare different Oq values by requiring the same 
initial gravitational potential $ (below the curvature scale). For scale invariant n = 1 A 
models, this normalization is equivalent to that determined by large scale anisotropics, since 
the SW effect dominates all but the lowest multipoles. Therefore the approximate relative 
amplitude of matter fluctuations can be directly read off from this plot. For open models, 
this is not true due to a more significant ISW effect and curvature effects at large scales 
which relate the potential to the initial power spectrum. 

Figure 17. The isocurvature matter transfer function (fJo = 0.2, h = 0.5, numerical). 
The baryon perturbations will have a prominent peak at the maximal Jeans scale since 
perturbations grow as D(a) outside this scale, with a fc^ ~3K tail, and are suppressed inside 
it. The acoustic oscillations damp away in the highly ionized case since last scattering is 
delayed and the diffusion length grows. This leaves a flat small scale tail in the transfer 
function. Note also that in the low ionization scenarios, the Jeans length may not grown 
to its maximum matter dominated value by last scattering leading to a smaller scale for 
the peak. The growth suppression due to A is less significant than that from curvature 
domination. 

Figure 18. The full adiabatic photon power spectrum. The logarithmic contributions to 
the anisotropy in £ and k [see equation (9)] are plotted here. Whereas in the Jlo = 1 case 
only the SW ridge and acoustic peaks are prominent {top left and close up, top right), the A 
and open cases show more complicated structure due to the ISW effect. Depending on the 
initial weightings, represented here as |CyiP oc k, certain features may be emphasized over 
others. Notice the A ISW effect at low £ and intermediate k and the comparatively small 
open SW contribution at the foot of the ISW ridge. 
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Figure 19. Large angle anisotropy dependence on the initial power spectrum jC/p ex 
fc" or |CaP oc fc" for isocurvature and adiabatic scenarios respectively, (fio = 0-lj h = 
0.5, standard recombination). Notice that for red spectra, geometric and/or cosmological 
constant effects play a role in determining the anisotropy whereas for very blue spectra, 
i'^Ci ex £'^ for all models. Isocurvature models with n ~ — 1 to fit large scale structure will 
thus not be extremely sensitive to open or A dominated universe effects. The normalization 
here is arbitrarily set at the quadrupole. 

Figure 20. The full open isocurvature photon power spectrum for |C/p ex A:", (a) Cur- 
vature scale weighted n — — 3. The lack of a Poisson cut off in the isocurvature potential 
makes the nature of the open universe eigenfunctions apparent. The projection ridge crosses 
minimum eigenvalue k = \/ —K (front edge) at roughly £ ~ 10 corresponding to the fact 
that the lowest eigenfunction still has structure only around the curvature scale. This leads 
to the cutoff to low multipolcs depicted in Fig. 19. (b) Small scale weighted n = —1. The 
main projection ridge does not dominate the anisotropy at the low order multipoles. Power 
from smaller physical scales (high k) bleeds in to boost the anisotropy. Thus anisotropics 
do not rise as rapidly with £ as predicted from the one to one conversion of k onto i. For 
this model ricff ~ 2 at large angles. 

Figure 21. Intermediate to small scale anisotropics, (a) Adiabatic models. Projection 
of the sound horizon at last scattering onto sky today determines the angular scale of the 
"Doppler peaks" (c./. Fig. 11). The sound horizon is the same physical scale for open and 
A models with fixed fig but geodesic deviation makes it correspond to a smaller angle in 
the open case. Compared with the flat case, the A model has a somewhat smaller angular 
scale due to the imperfect cancellation between the increase in the age of the universe today 
and at last scattering, (b) Isocurvature models. Anisotropics in the standard recombination 
scenario {xe ~ if ;z < 1000) produce far to large fluctuations on the arcminute scale due to 
the steeply rising spectrum. Reionized models have their adiabatic fluctuations damped out 
by photon diffusion and a cancellation suppressed Doppler effect. Notice that large angle 
anisotropics are immune to ionization history effects for the open case but not for the A 
case. This and the difference in the damping scale is mostly due to the projection effect. 

Figure 22. Isocurvature temperature power spectrum. In this fully ionized x^ = 1, low 
flo = 0.1 h = 0.5 model, the ISW effect makes a contribution equal to and with the 
same scale dependence as the cancelled Doppler (plus SW) term (DSW). The second order 
Vishniac term (V) dominates at small scales. The analytic approximation (solid) fails at 
large scales where cancellation arguments are not applicable. 

Figure 23. Radial eigenfunctions of an open universe X^ix)- (a) The isotropic £ = 
function for several values of the wavenumber ly = kj \f--K. The zero crossing moves out to 
arbitrarily large scales as :/ ^- 0, reflecting completeness. However, even as this "effective 
wavelength" becomes infinite, the function retains prominent structure only near the cur- 
vature scale X- ^ random superposition of these low v modes cannot produce more than 
exponentially decaying structure larger than the curvature scale, (b) Low order multipoles 
in the asymptotic limit j/ ^ 0. If most power lies on the curvature scale, the ^-mode corre- 
sponding to the angle that the curvature radius subtends will dominate the anisotropy. The 
normalization is appropriate for comparing contributions to the anisotropy £(2i?+ l)C^/47r. 
Also shown is the location of the horizon x = VoV—K for several values of Dq- If contribu- 
tions to the anisotropy come from a sufficiently early epoch, the dominant ^-modc will peak 
at this value. 
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